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Abstract 

Multiple orthogonality is considered in the realm of a Gauss-Borel factorization problem for a semi-infinite 
moment matrix. Perfect combinations of weights and a finite Borcl measure are constructed in terms of M-Nikishin 
systems. These perfect combinations ensure that the problem of mixed multiple orthogonality has a unique solution, 
that can be obtained from the solution of a Gauss-Borel factorization problem for a semi-infinite matrix, which 
plays the role of a moment matrix. This leads to sequences of multiple orthogonal polynomials, their duals and 
second kind functions. It also gives the corresponding linear forms that are bi-orthogonal to the dual linear forms. 
Expressions for these objects in terms of determinants from the moment matrix are given, recursion relations are 
found, which imply a multi-diagonal Jacobi type matrix with snake shape, and results like the ABC theorem or 
the Christoffel-Darboux formula are re-derived in this context (using the factorization problem and the generalized 
Hankel symmetry of the moment matrix). The connection between this description of multiple orthogonality 
and the multi-component 2D Toda hierarchy, which can be also understood and studied through a Gauss-Borel 
factorization problem, is discussed. Deformations of the weights, natural for M-Nikishin systems, are considered 
and the correspondence with solutions to the integrable hierarchy, represented as a collection of Lax equations, 
is explored. Corresponding Lax and ZaJtharov-Shabat matrices as well as wave functions and their adjoints are 
determined. The construction of discrete flows is discussed in terms of Miwa transformations which involve Darboux 
transformations for the rrmltiple orthogonality conditions. The bilinear equations are derived and the T-function 
representation of the multiple orthogonality is given. 
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1 Introduction 

The topic of multiple orthogonality of polynomials is very close to that of simultaneous rational approximation (si- 
multaneous Pade aproximants) of systems of Cauchy transforms of measures. The history of simultaneous rational 
approximation starts in 1873 with the well known article [22J in which Ch. Hermite proved the transcendence of the 
Euler number e. Later, around the years 1934-35, K. Mahler delivered at the University of Groningen several lectures 
[?7] where he settled down the foundations of this theory. Meanwhile, two of Malher's students, J. Coates and H. Jager, 
made important contributions in this respect (see (I3j and j23| ). In the case of Cauchy transforms, the simultaneous 
rational approximation definition may be written in terms of multiple orthogonality of polynomials as follows. Given 
an interval A C M of the real line, let A^(A) denote all the finite Borel measures which have support, supp(-) with 
infinitely many points in A, where they do not change sign. Fix € A^(A), and let us consider a system of weights 
w = {wi, . . . , Wp) on A, with p € N. (In this paper a"weight" on an interval A is meant to be a real integrable function 
defined on A which does not change its sign on A.) Fix a multi-index i? ~ {vi, . . . , Vp) G Z^, Z_|- = {0, 1,2,.. .}, and 
denote = vi + ■ ■ ■ + Vp. There exist polynomials, Ai, . . . , Ap, not all identically equal to zero which satisfy the 
following orthogonality relations 

r P 

x'Y.Aa{x)wa{x)dfi{x)^0, degA,<i.,-l, j = 0,...,|i7|-2. (1) 

•^A a=l 

Analogously, there exists a polynomial B not identically equal to zero, such that 

[ x'B{x)wb{x)dn{x)^0, degB<|i7|, j = 0, . . . , z^h - 1, 6=1,..., p. (2) 

J A 
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The resulting polynomials are said to be of type I and type II, respectively, with respect to the combination (/i, w, u) 
of the measure /i, the systems of weights w and the multi-index v. When p ~ 1 both definitions coincide with that of 
the standard orthogonal polynomials on the real line. The existence of a system of polynomials {Ai, . . . ,Ap) and a 
polynomial B defined from ([T]) and ([2| respectively, are ensured because in both cases finding the coefficients of the 
polynomials is equivalent to solving a system of linear homogeneous equations with + 1 unknown coefficients. 
From the theory of orthogonal polynomials we know that when p = 1 each polynomial Ai = B has exactly degree | iy\ = 
vi \ unfortunately if p > 1 that is not true in general. For instance, let us take a system of weights w = (wi, wi, . . . , wi), 
in this case the solution vector space has dimension bigger than one, and we can find two solutions which are linearly 
independent. Hence, there is at least an a € {1, . . . .p} such that degAa < Va — \ and deg_B < \v\. Given a measure 
/i G A1(A) and a system of weights w on A a multi-index v is called type I or type II normal if deg must equal to 

— 1, a = 1, . . . ,p, or deg_B must equal to \v\ — 1, respectively. When for a pair (/i, w) all the multi-indices are type 
I or type II normal, then the pair is called type I perfect or type II perfect, respectively. The concepts of normality 
and perfectness were introduced by Malher (see Malher's, Coates' and Jager's articles cited above). 

Multiple orthogonal of polynomials have been employed in several proofs of irrationality of numbers. For example, 
in fTO|, F. Beukers shows that Apery's proof (see [8^) of the irrationality of ^(3) can be placed in the context of 
a combination of type I and type II multiple orthogonality, which is called mixed type multiple orthogonality of 
polynomials. More recently, mixed type approximation has appeared in random matrix and non-intersecting Brownian 
motion theories (see, for example, [TT], [H] and |"25J). A formalization of this kind of orthogonality was initiated by 
V. N. Sorokin (36) . He studied a simultaneous rational approximation construction which is closely connected with 
multiple orthogonal polynomials of mixed type. Surprisingly, in |21j a Riemann-Hilbert problem was found for the 
theory of orthogonal polynomials, and later [3^ this result was largely extended to type I and II multiple orthogonality. 
In [141 mixed type multiple orthogonality was analyzed from this perspective. 

In order to introduce multiple orthogonal polynomials of mixed type we need two systems of weights wi = 
(wi_i, . . . , wi^pj) and W2 — iw2,i, ■ ■ ■ ,W2,p2) where pi,P2 G P^, (as we said a set of functions which do not change 
their sign in A), and two multi-indices vi = (i^i^i, . . . , vi^pi) G It^^ and V2 = (i^2,i, • • • , '^2,p2) ^ ^+ with — \v2\-\' 1- 
There exist polynomials Ai, . . . , , not all identically zero, such that deg^s < '^i,s which satisfy the following 
relations 

I- pi 

I ^S^Aa{x)wi,a{x)w2.b{x)x^d^i{x) ^Q, j = 0, . . . , Z/2,6 - 1, 5=1,...,P2- (3) 

They are called mixed multiple-orthogonal polynomials with respect to the combination {fi^wi,'W2, 1^1,1^2) of the 
measure n, the systems of weights wi and W2 and the multi-indices vi and i'2- It is easy to show that finding the 
polynomials Ai, . . . ,Ap-^ is equivalent to solving a system of \iy2\ homogeneous linear equations for the \vi\ unknown 
coefficients of the polynomials. Since — |i'2| + 1 the system always has a nontrivial solution. The matrix of 
this system of equations is the so called moment matrix, and the study of its Gauss-Borel factorization will be the 
cornerstone of this paper. Observe that when pi = 1 we are in the type II case and if p2 = 1 in type I case. Hence 
in general we can find a solution of ([3| where there is an a e {1, . . . such that deg Aa < i^im — ^- When given a 
combination (//, i/Ji, 7112) of a measure fj, G A^(A) and systems of weights wi and W2 on A if for each pair of multi- 
indices {vi, V2) the conditions (|3| determine that deg = vi a — 1, a = 1, . . . then we say that the combination 
{fi,wi,'W2) is perfect. The concept of perfectness will be rigorously introduced in Definition [2] 

The seminal paper of M. Sato and further developments performed by the Kyoto school through the use of 
the bilinear equation and the r-function formalism [l6]-[18], settled the basis for the Lie group theoretical description 
of integrable hierarchies, in this direction we have the relevant contribution by M. Mulase [30] in which the factor- 
ization problems, dressing procedure, and linear systems were the key for integrability. In this dressing setting the 
multicomponent integrable hierarchies of Toda type were analyzed in depth by K. Ueno and T. Takasaki [37| . See 
also the papers [H] and [51] on the multi-component KP hierarchy and [2H] on the multi-component Toda lattice 
hierarchy. In a series of papers M. Adler and P. van Moerbeke showed how the Gauss-Borel factorization problem 
appears in the theory of the 2D Toda hierarchy and what they called the discrete KP hierarchy [l]-|5]. In these 
papers it becomes clear -from a group-theoretical setup- why standard orthogonality of polynomials and integrability 
of nonlinear equations of Toda type where so close. In fact, the Gauss-Borel factorization of the moment matrix 
may be understood as the Gauss-Borel factorization of the initial condition for the integrable hierarchy. To see the 
connection between the work of Mulase and that of Adler and van Moerbeke see [H] . Later on, in the recent paper 
[6] , it is shown that the multiple orthogonal construction described in previous paragraphs was linked with the multi- 
component KP hierarchy. In fact, for a given set of weights (wi,?«2) and degrees {1^1, i'2) the authors constructed a 
finite matrix that plays the role of the moment matrix and, using the Riemann-Hilbert problem of |14| , where able to 
show that determinants constructed from the moment matrix were r-functions solving the bilinear equation for the 
multi-component KP hierarchy. However, there is no mention in that paper to any Gauss-Borel factorization in spite 
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of being the multicomponent integrable hierarchies connected with different factorization problems of these type. For 
further developments on the Gauss-Borel factorization and mult i- component 2D Toda hierarchy see [7 and 29J. 

This motivated our initial research in relation with this paper; i.e., the construction of an appropriate Gauss- 
Borel factorization in the group of semi-infinite matrices leading to multiple orthogonality and integrability in a 
simultaneous manner. The main advantage of this approach lies in the application of different techniques based on the 
factorization problem used frequently in the theory of integrable systems. The key finding of this paper is, therefore, 
the characterization of a semi-infinite moment matrix whose Gauss-Borel factorization leads directly to multiple 
orthogonality. This makes sense when factorization can be performed, which is the case for perfect combinations 
(/I, iSi, z«2), which allows us to consider some sets of multiple orthogonal polynomials (called ladders) very much 
in the same manner as in the (non multiple) orthogonal polynomial setting. The Gauss-Borel factorization of this 
moment matrix leads, when one takes into account the Hankel type symmetry of the moment matrix, to results like: 
1. Recursion relations, 2. ABC theorems and 3. Christoffel-Darboux formulas. The first two are new results while 
the third is not new, as it was derived from the Riemann-Hilbert problem in |14) . However, our derivation of the 
Christoffel-Darboux formula is based exclusively on the Gauss-Borel factorization, and its uniqueness and existence 
for the multiple orthogonality problem are the only requirements. Thus, it is sufficient to have a perfect combination 
{fj,,Wi,'W2), and there are examples of this type which do not have a well defined Riemann-Hilbert problem in the 
spirit of [Ti] . 

When we seek for the appropriate integrable hierarchy linked with multiple orthogonality we are lead to the 
multicomponent 2D Toda lattice hierarchy which extends the construction of the multicomponent KP hierarchy 
considered by M. J. Bergvelt and A. P. E. ten Kroode in [9]; not to the multicomponent 2D Toda lattice hierarchy as 
described in '37 or |28| . In the spirit of this last mentioned articles, and complementing the continuous flows of the 
integrable hierarchy, we also introduce discrete flows, that could be viewed as Darboux transformations, and which 
correspond to Miwa transformations implying the addition of a zero/pole to the set of weights. Moreover, the Hankel 
type symmetry is related to an invariance under a number of flows, and to string equations. Bilinear equations can be 
derived from the Gauss-Borel factorization problem and moreover the r-function representation is available leading 
to a bridge to the results of 6J in which no semi-infinite matrix or Gauss-Borel factorization was used. 

This paper is divided into three sections, §1 is this introduction which contains §I.I| in where we review the 
application of the LU factorization of the moment matrix to the theory of orthogonal polynomials in the real line. 
Next, §2 is devoted to the presentation of the moment matrix and the discussion of the Gauss-Borel factorization. In 
this form we obtain perfect systems in terms of Nikishin systems, determinantal expressions for the multiple orthogonal 
polynomials, their duals and second type functions, bi-orthogonality for the associated linear forms, recursion relations, 
ABC type theorems and the Christoffel-Darboux formula. Flows and the integrable hierarchy are studied in §3 in 
which an integrable hierarchy a la Bergvelt-ten Kroode is linked with the multiple orthogonality problem. We not 
only derive from the Gauss-Borel factorization the Lax and Zakarov-Shabat equations, but also we introduce discrete 
integrable flows, described by Miwa shifts, or Darboux transformations, and also construct an appropriate bilinear 
equation. Finally, we find the r functions corresponding to the multiple orthogonality and link them to those of f6j. 
At the end of the paper, we have added two appendices: the first one collects the more technical proofs of the results 
in this paper. In Appendix [B] we consider discrete fiows for the case of a measure ^ with unbounded support supp/x. 

1.1 The Gauss-Borel factorization of the moment matrix and orthogonal polynomials 

Here we discuss how the LU factorization of the standard moment matrix g = (Jx'+M/i) of a constant sign finite 
Borel measure fj, leads to traditional results in the theory of orthogonal polynomials, namely recursion relation and 
Christoffel-Darboux formula. In spite that these results are well established we repeat them here because in their 
derivation is encoded the set of arguments we will use in the multiple orthogonality setting. In the forthcoming 
exposition it will become clear the LU factorization approach is just a compact way of using the orthogonality 



relations. 

The moment matrix can be written as the following Grammian matrix 

9 = X{x)x{x)'^ d^i{x) 



in terms of the monomial string xi^) '■— (1; 2;, x^. 
The Borel-Gauss factorization of g is 
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The reader should notice that 



• It makes sense whenever the truncated moment matrix gl'l — {gi.j)Q<i.j<:i is an invertible matrix for any I = 
1,2,.... If the factorization exists it is unique. 

• Although the truncated matrices g^'l are invertible it can be shown that g itself is not invertible. 

• The matrix product of S^^ with S involves only finite sums, but if we reverse the order of the factors we get 
series (with an infinite number of summands). 

Given the factors S and S we consider the following polynomial strings, the semi-infinite vectors, 
P:^Sx^ (Po, Pi,--. P := {S^Yx = (^o, Pi,--- V- 

The families of polynomials {Pi}'^q and {Pk}'kLa are biorthogonal: 

P{x)P{xyd^i{x) = / Sxix)xixV S-'dflix) j x{x)x{Y d^i{x)S~^ 

Pi{x)Pk{x)dfi{x) = Sisk- 



in this simple proof relics the basic connection between orthogonality and the LU factorization, which we consider as 
the very same thing dressed in different manners. From the above orthogonality we conclude that 



Pi{x)x'dfi{x)^0, i = 0,...,/-l, 
Pi{x)x^ dfiix)^0, j=0,...,Z-l, 



(4) 



and we also have that Pi{x) and Pi are Z-th degree polynomials where Pi is monic and Pi satisfies J x^ Pi{x)d^{x) = 1, 
i.e. we have type II and type I normalizations. Given that the moment matrix is symmetric, g = and the 
uniqueness of the LU factorization we deduce that S = H{S^^)^ , with H = diag(/io, hi,. . . ); i.e., g = S~'^H{S~^)^ 
and the factorization is a Cholesky factorization (but this does not extend to the multiple orthogonal case) . Therefore 
Pi = h^^Pi so that J Pi{x)Pk{x)dii{x) = hiSi,k, and {i^ij^o ^ family of monic orthogonal polynomials with respect 
to the measure fj.. 

Considering the orthogonality relations as a linear system for the coefficients of the polynomials one concludes that 
polynomials and their duals can be expressed as 



= x^'' - 




91,1 ■ 


■ gi.i-i) 






- Su (0 





- - 


1) (g['+ii) 


-i^[^+i] 








/ 50,0 


50,1 


5o,;-i 


1 


1 

^ det ^['l 


det 


51,0 


51,1 


5i,i-i 


X 






5/-i,o 


.9/-i,i 


gi-i,i-i 






\ 9i,o 


.9i,i 


gi.i~i 


x' 



> 1, 



/ 

and similar expressions for the dual polynomials. We are now ready to get the recursion relations for orthogonal 
polynomials: 

• First, we notice that the moment matrix 5 is a Hankel matrix, gi+ij = g^.j+i, which in terms of the shift matrix 
'0 1 ... 



A 



0010 
0001 



can be written as Kg = gbJ 
• Second, we observe the eigen- value property Ax(a;) ~ xxix). 

^From this symmetry property it follows, by contradiction, that the moment matrix is not invertible; i.e. the assumption of the existence 
of leads to g~^A = A^g"^, and therefore the first row and column of are identically zero, so that is not invertible. 
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• Third, we introduce the LU factorization to get AS^^S S^^SA^ ^ SAS^'^ = S'A^S'^^ -. J. From this last 

(ao 1 ... \ 
6i ai 1 ... \ 
62 02 1 is a tridiagonal matrix, i.e. a Jacobi matrix. 

• Finally, we notice that the polynomial strings are eigenvectors of the Jacobi matrix: JP{x) = SAS~^ Sx{x) = 
SAx{x) = Sxx{x) = xP{x)\ i.e., the recursion relations xPfe(x) = Pk+i{x) + a^Pkix) + bkPk~i{x), k > 0, hold. 

We now consider the Aitken-Berg-CoUar (ABC) theorem (here we follow the nomenclature used [55]) foi' orthogonal 
polynomials. First we introduce the Christoffel-Darboux kernel and therefore we consider 

■^['1 =M{0...,a;'-i}, n = { J2 cix^cieRy (-hW)"^ = { ^ QP'(x),QeM} 

0<i<oo ;<fc<oo 

and the resolution of the identity H = "H''' ® ('HI'I)^, with the corresponding orthogonal projector tt*^'' such that 
ker TT^') = ('H^)-'- and Ran tt^'^ — "HI'I. Then, the Christoffel-Darboux is defined as 

i-i i-i 
xW(x,y) :=J2Pkiy)Pk{x)=Y.h-'Pk{y)Pk{x), 

k=0 k=0 

which, according to the bi-orthogonality property, gives the following integral representation of the projection operator 

^('V)(y) = / K^Hx,y)f{x)d^i{x), \ff e n, 



Any semi-infinite vector v can be written in block form as follows 

is the finite vector formed with the first I coefficients of v and the semi-infinite vector formed with the 

remaining coefficients. This decomposition induces the following block structure for any semi- infinite matrix. 



gll] 


gll,>l] 


gl>l,l] 





Given a factorizable moment matrix g we have 

gll] ^ ^ (^W)-i^ (5-i)[^'l = {S^^'^)-\ 

The Christoffel-Darboux kernel is related to the moment matrix in the following way 

which is a consequence of the following identities 

i^['l(x,y) = (nWp(a;))^(nWp(y)) 
= X^(a:)5-inW5x(2/) 

= (x) (ni'i ^-^nw ) (ni'i sn^ )x(y) 

= ix^Hx)f{gVy)-\W^yy 

The relations 

(.9W)-1aW - (AW)T(gW)-i = (gW)-i(5['^>'l(A['^^'l)T - Al^'^'ly^'^'l) (.gW)"! 
follow from the block equation 

j^Wgli] ^ j^V,>i]g[>i,i] = g['l(AW)^ +5l''^'l(A[''^'l)^. 
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We also have 

where {e^j^Q is the canonical linear basis of T-L. With all these at hand we deduce 
so that, 

That using the determinantal expressions for the polynomials presented before leads to the Christoffel-Darboux formula 

(x- y)xW(x,y) = hi\{Pi{x)Pi^^{y) - Pi^^{x)Pi{ij)). 

2 Multiple orthogonal polynomials and Gauss— Borel factorization 
2.1 The moment matrix 

In this section we define the moment matrix in terms of a measure /i G A^(A) and two systems of weights wi 
and W2 on A C M, as well as corresponding compositions (the order matters) fii = (ni^i, . . . , ni^p^^) € and 
1 ^2,p2) G f^^^ l38J. We will consider multi- indices of positive integers n — (tt-i, . . . , rip), where p G N and 
ria e Z+, a = 1, . . . ,p and define |n| := rii + • • • + Up. Following (21 [35] we observe that any z e Z+ := {0, 1,2,...} 
determines unique non-negative integers q{i),a{i),r{i), such that the composition 

i ^ q{i)\n\ + ni ^ h na(i)-i + ''l*), < r(i) < n<j(j), (5) 

holds. Hence, given i there is a unique k{i) with 

k{i) ^ q{i)na{i)+r{i), < r(i) < n<j(j)- (6) 

Let us introduce the function integer part function [•] : — > Z-)_, [a;] max{j/ e "^+,1/ < x}. Combining ([s]) and ^ 
we can obtain a formula which expresses explicitly the dependence between the quantities i, k and a 



|"|-«a)+"-iH hria-i + fc- (7) 



Let M°° denote the vector space of all sequences with elements in M. An clement A G may be interpreted as a 
column semi-infinity vector as follows 

A = (A(o),A(i),...)^, A«eE, J =0,1,.... 

We consider the set {ej}j>a C M°° with 

ej = (OX^^, 1,0,0,...)^. 

Here (•)^ denotes the transposition function on vectors and matrices. Analogously, we denote by the set 

of all sequences of vectors with p components and observe that each sequence which belongs to can also be 

understood as semi-infinity column vector: given the vector sequence {vo,vi, . . . ) with vj = . . . j'l'j.p)^ we have 
the corresponding sequence in R°° given by (vo,i • ■ • , vo,p, . • . , wi,p, . . . ); i.e., R°° = {Rp)°° . Therefore, we consider 
also the set {ea{k)}a=i,....p C where for each pair (a, fc) € {!,..., p} x Z+ ea{k) — ei(^ka) and the function 

k=0,l,'... 

i{a,k) £ Z+ satisfies the equality ([t]). 

Now, we are ready to introduce the monomial strings 

X„:=f:e.(fc)z^ Xi'^ = |f'^' x:: = ^-V(^-^). (8) 
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These vectors may be understood as sequences of monomials according to the composition n introduced previously. 
We also define the following weighted monomial string 



C:=j2xaw^, C^'^ = w,ii)z'^'\ (9) 

a=l 

which is a sequence of weighted monomials for each given composition ft. Sometimes, when we what to stress the 
dependence in the composition we write Xn,a,Xn and ^fj. Given the weighted monomials ^ft^ and ^ft^, associated to 
the compositions Hi and fi2, we introduce the moment matrix in the following manner 

Definition 1. The moment matrix is given by 

' ^nAx)^nA^)'dtx(x). (10) 



In terms of the canonical basis {-E/.j} of the linear space of scmi-infinitc matrices and for each pair G 
we consider the binary permutations or transpositions Wij = Eij + Ej^i. Observe that n^j = I and therefore n^J = 
TTij. Given two transpositions Wij and TTfc,; the permutation endomorphism corresponding to its product is well 
defined TrijiTk.i = '^k.i'^i.j- Taking a sequence of pairs / — {{is, js)}sez+, is,js G 2+, we introduce the permutation 
endomorphism tt/ as the infinite product ttj = nsez+ '^iaJs' with ttjttJ = ttJttj = I. Given two compositions, n',n, 
there exists a permutation tt^j/^^j such that ^^5/ = nff^ft^ft through a permutation semi- infinite matrix as just described. 

The change in the compositions is modeled as follows 

Proposition 1. Given two set of weights we = {w(^i, ... jWe^p^) and compositions ne and n'^, I = 1,2, there exist 
permutation matrices tt^^^s^ such that 

Proof. For any set of weights w = wi, . . . ,Wp and two compositions ft and n' we have that the corresponding vectors 
of weighted monomials are connected, 

trough a permutation somi-infinito matrix; i.e, ttT, - = Tf^^/^-j- Therefore, the announced result follows. □ 

For the sake of notation simplicity and when the context is clear enough we will drop the subindex indicating the 
two compositions and just write g for the moment matrix. Let us discuss in more detail the block Hankel structure 
of the moment matrix. For each pair G there exists a unique combination of three others pairs {qi,q2) G Z^, 
(ai,a2) G {1, . . . ,pi} X {1, . . . ,p2} and (ri,r2) G {0, . . . ,ni,ai - 1} x {0, . . . ,n2,02 - 1}> such that 

^ = al|«l|+?^l,l + •••+?^l,al-l+n and j = g2|ri2| + ^2,1 + • • • + n2,a2-i + ^2- 

Hence taking k(, = qene^ae +re, i = 1,2, the coefficients gij G IR of the moment matrix g = (gij) have the following 
explicit form 



J x'''+''-Wi,aAx)W2,aAx)dn{x). (12) 



Observe that pairs {ki,ai) and (^2,02) are univocally determined by i and j respectively. 

Before we continue with the study of this moment matrix it is necessary to introduce some auxiliary objects 
associated with the vector space First, we have the unity matrix I — J^kLo^^eJ and the shift matrix A := 

X^^Lo efceLi- We also define the projections Ill'l := X^tlnefee^, and with the help of the set {ea{k)}a=i,...,p we 

/c=0,l,... 

construct the projections Ila := Y^'^^o ^a{k)ea{k)''' with J2a=i ~ ^' ^^'^ 

Pi :=diag(I„i,0„2,...,0„p), P2 := diag(0„i,I„2, . . . , 0„^), ... Pp := diag(0„i,0„2, . . . ,I„^), (13) 
where I„_^ is the rig x rig identity matrix. Finally we introduce the notation 

x^ :=a;"iPi + ---+a;"''Pp = diag(x"iI„,,...,a;"''I„J :M^mI"IxI"I. (14) 
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For a better insight of the moment matrix let us introduce the following ni_a x f^2,b matrices 



/ 



ma,b{x) = Wi^aix)w2,bix) 



1 

X 



X 



rai,<i+n2,!, — 2 



,Pl, 
,P2, 



(15) 



in terms of which we build up the following \ni \ x |ri2|-matrix 

/ mi,i mi, 2 • • • ™i,p2 \ 



7712,1 "l2,2 



™2,p2 



\r77pi,i 



p|ni|x|n2| 



(16) 



1pi,2 



Then, the moment matrix g has the following block structure 
9 '■= {Gij)i,j>o G Gij 



(17) 



Fix now a number I G N and consider the pair (1,1 + 1). There exists a unique combination of pairs (gi, 52) G 
(ai,a2) e {1, • • • X {1, . . . ,p2} and (ri,r2) G {0, . . . ,77i,q^ - 1} x {0, . . . , n2,a2 - 1}, such that 



+ ^^l.al-l + T"! and ^ + 1 = (72|"2| + "2,1 + 



"2,02 



^2. 



Given the compositions tii and 772 we introduce the degree multi-indices i^i G and 1^2 G [9] where for each 
: 1, 2, we have 



which satisfy 

kfXi + l) = vi^a,_(i){i), = 
and consider the I x (Z + 1) block matrix F; from g 



1, 



l7f (« + l\nt\) = 77£(7) + /77£, 



5i,o 



5o,i 
51,1 



5o,i \ 



(18) 



(19) 



(20) 



Let us study the homogeneous system F;a;/_|_i = Oi, where xi+i G and Oi+i is the null vector in M}^'^. Taking into 
account Fj's structure (12 1, we see that such equation is exactly the expression of the orthogonality relations (|3|. We 
can see now that for each I G N the existence of a system of mixed multiple orthogonal polynomials [A'^\ . . . ,Ap}^ 
is ensured; that is because F; in (20l is a Z x (Z + 1) matrix, and the homogeneous matrix equation Fjccj+i = 0/, 
which is satisfied by the coefficients of the polynomials corresponds to a system of I homogeneous linear equations 
for I + 1 unknown coefficients. Thus, the system always has a non-trivial solution. Obviously, (^i'\ . . . is not 

univocally determined by the matrix equation F/ccj+i — 0/ or equivalently by the orthogonality relations ([s]), because 
its solution space has at least dimension 1. Hence, the appropriate question to consider is the uniqueness question 
without counting constant factors, or equivalently if the solution space has exactly dimension 1. In terms of F; the 
question becomes: Does F; have rank 17 In order to have a positive answer it is sufficient to ensure that the I x I 
square matrix 



/ 9o,o 
91,0 



9oA 
91,1 



\5;-i,o 9i-i,i 



90.1-1 \ 
91,1-1 

91-1,1-1/ 



I > 1, 



(21) 



is invertible, where t/W results from F; after removing its last column. It is easy to prove that such condition is 
equivalent to require that all possible solutions of ^ satisfy deg Ap.^ — vi^p^ — 1. Obviously this requirement is 



ensured when the polynomials (A^'-*, . . . , Apl) fulfill degAj 



1, J = 1, 
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2.2 Perfect combinations and Nikishin systems 

We introduce the concept of a perfect combination. 

Definition 2. A combination {fi,wi,'W2) of a measure fj, G A^(A) and two systems of weights wi and W2 on A C M 
is said to be perfect if for each pair of multi-indices {vi,V2), with \vi\ = \v2\ + 1 the orthogonality relations ([3| imply 
that deg^d — vi^a ~ ol = 1- • • ■ jPi- 

For a perfect combination wi,^'W2) and any given I G the solution space of the equation TiXi^i — 0/ is one- 
dimensional. Then, we can determine a unique system of mixed type orthogonal polynomials (Ai, . . . , Ap^) satisfying 
([3]) requiring for oi G {!,... pi\ that Aa-^ monic. Following 14] we say that we have a type II normalization and 
denote the corresponding system of polynomials by A^a^'°'^\ j — 1, . . . Alternatively, we can proceed as follows, 
since the system of weights is perfect from ([3| we deduce that 

pi 

C^^-O^y^ Aa{x)wi^aix)w2.b2{x)dKx) + 0. 



Then, we can determine a unique system of mixed type of multi-orthogonal polynomials (^^^'''^'', . . . , a!^^^^^) imposing 
that 

' A^l-''^\x)wx^a{x)W2.hWil[x) = 1, 



Pi 

,1(1,62 

a=l 



which is a type I normalization. We will use the notation aIV'I^I and ylf!.'''-\ to denote these multiple orthogonal 
polynomials with type II and I normalizations, respectively. 

A known illustration of perfect combinations (/i, wi, W2) can be constructed with an arbitrary positive finite Borel 
measure /i and systems of weights formed with exponentials: 

(e^^^...,e>-), 7*^7,, » ^ J, *,j = l,...,p, (22) 

or by binomial functions 

((l-z)"S...,(l-z)"-), a,-a, ^Z, z^j, *,j = l,...,p. (23) 

or combining both classes, see {5T1. Recently a wide class of systems of weights where proven to be perfect j20] : 
these systems of functions, now called Nikishin systems, were introduced by E.M. Nikishin |31j and initially named 
MT-systems (after Markov and Tchebycheff ). 

o 

Given a closed interval A let A be the interior set of A. Let us take two intervals A^ and A^ whose interior sets 

o o 

are disjoint, i.e. A^ n A,3 = 0. Set two measures \Xa G M.{^a) and [ip G A^(A^) such that the measure (/Iq,, with 
the following differential form 

d{fla,tJ'P){x) = J dfia(x) = flp{x)dHa{x), 

is a finite measure, that implies that (/ic, fip) G A4{Ao,). The function p,p denotes the Cauchy transform corresponding 

o o 

to fip. Let us consider then a system of p intervals Ai, . . . , Ap such that A^ n A^+i = 0, j G {1, ... ,p — 1}. Take p 
measures fij G A4{Aj), which for each j = 1, . . . ,p — 1, the measure (/ij, cTj+i) belongs to A4{Aj). So the system of 
measures (^o, ■ • ■ : Cp) where 

Cl==Ml, C2=(Ail,M2), C3 = (Ai2,Ai3)) = (Ail,Ai2,/^3), Cp = • ■ • ,Aip), 

is the Nikishin system of measures generated by the system (cti, . . . , am)- So we denote (Ci, . . . , Cp) ^ Af (ui, . . . , Cp) . 

Actually, in [2^ the authors shown perfectness for combinations of Nikishin systems where intervals Ai, . . . , Ap are 
bounded and for each j G {1, . . . ,p — 1} the intervals Aj and Aj^i are disjoint. The same authors have communicated 
to us that they were able to prove a generalization of this result to unbounded intervals such that A^ n Aj+i 7^ 0. 
Consequently, in what follows we assume such generalization. 
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As we have seen, general Nikishin systems have an intricate structure; therefore, in order to make easy the reader 
we focus on a "simple" class of Nikishin systems which we call M-Nikishin systems. Set the interval Ai ~ [0, 1] and 
let 7Mo(Ai) C A^(Ai) denote the set of measures in A^(Ai) such that if cr e A^o(^i) then, the function 

\Aa{t)\ 



a(x) := 



satisnes lim 




= lim / 


\ — tx 

o 







1-tx 



< +00, where Ai = (0,1) 



(24) 



The constraint (24) guarantees that the function a is a weight in compact intervals in {—00, 1]. As (1 — tx) does not 
vanish for {t, a;) e Ai x (C\ [1, +00)) we deduce that 1/(1 — ta:) is a continuous function in a; for t e Ai. Therefore, we 
conclude that 5 is a holomorphic function on C \ (1, +00), having a continuation as continuous function in 1. Taking 
into account that a does not vanish in C \ (1, +00) and that it takes real values on M \ (1, +00) = (—00, 1], we deduce 
that it is a continuous weight on (—00, 1]. Observe that 



(7{x) = 



da{t) 
1 - tx 



Cda{\/0 



[l,+oo) 



MO 

x~C 



(25) 



[l,+oo) ^ C 

is the Cauchy transform of another measure n G A1([l, +00)), such that |/i(l)| = |ct(1)| < +00. 

Given two measures CTq G A^o(Ai),ct^ G Alo(^i) we define a third one as follows (using the differential notation) 

da;3(C) 



d[CTQ,cr/3](a;) = ap{x)daa{x), 



ap{x) = 



l-a;C 



As Up is a continuous weight on Ai we conclude that [ctq,, ap\ G A^o(^i)- If take a system of measures (ci, . . . , CTp) 
such that Gj G A^o(Ai)j J = ^, ■ ■ ■ ,P, we say that (si, . . . , Sp) = Jv[J\[{(Ji, . . . , (jp), where 

Sl=Cri, S2 = [0-l,Cr2], S3 = [CTI, [cr2,(T3]] = [o-l,CT2,0'3], ... Sp = [cri,0-2, ■ • • jCTp] (26) 

is the M-Nikishin system of measures generated by (cti, . . . , dp), with corresponding M-Nikishin system of fmictions 
given by w = (wi, . . . , Wp) (si, . . . ,Sp) = MU{cti, . . .,ap). 

Notice that s,; G Alo(^i) which implies that for each arbitrary compact subinterval of (— 00, 1] the system of 
functions w conforms a system of continuous weights. M-Nikihsin systems are included in the class of Nikishin 



systems. Taking into account the identity ([25| we see that the M-Nikishin system defined in (26 1 can be written as a 
classical Nikishin system. Let us take a system (fii, . . . , /Xp) where 

/ii = cri, dn2{x) = xda2{l/x), fJ.3 = as, fi2[p/2]-i ^ cr2[p/2]-i, Mlp/2]{x) ^ xa2[p/2]{l/x), 

and if p is odd fip — (Jp. Notice then 

Sl=Cl=Cri: S2 = C2 = (M1,M2), ■•■ Sp = Cp = (/^l7/^2, • ■ • 

Hence (si, . . . , Sp) = M.N {(Ti, . . . , dp) — A/'(^i, . . . , /ip) = (Ci, • . . , Cp)- Fixing two M-Nikishin systems of functions 
wi{x) = {si^i{x), . . .'si^p{x)) whose elements are weights on Aq = [—1, 1], and a measure fi G A4{Aq) we have at our 
disposal the perfect combination (/x, i/Ji, ^2)- We can also obtain a perfect combination [ijl,wi,W2) choosing wi and 
W2 between two different of the classes mentioned in ( 22 1 and ( 23 ) (not necessarily the same) . 

Proposition 2. The Taylor series at — corresponding to the functions Sj(C) and fj{C) ■— fogSj(C) converge 
uniformly to Jj and fj respectively on Ai, i.e. 



where Xij and tij 

Proof. For each j G {1, 
implies that 



Notice that 



i 

are constants. 



a: G Ac 



J = 1, ■ 



, ,p. 



(27) 



,p}, Sj is a holomorphic function on the open unitary disc centered on the origin. That 



Sjix) 



00 



i=0 



i=0 



lim 

1 

ce[o,i) 



1=0 



xG{|C|<l}, 
ds,(t) 



lim 



xeAi 



l-xt 



< +00, 



So the first equality in (27) is proved. The second one comes immediately from the fact that the functions sj do 
not vanish on Aq. That implies that X)i^o''*J^* bounded and therefore continuous. Hence we can proceed 

analogously as in the first equality. □ 
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2.2.1 The inverse problem 

Given the series 



Wj{x) 



i=0 



X e Ao, 



(28) 



we consider the problem of finding conditions over {Xij} such that the set of series {wj}^^-^ form a M-Nikishin system 
of functions. The reader should notice that Xi j = Si{tifi,ti^i, . . . ,tij) where Si is the i-th elementary Schur polyno- 
mial. Elementary Schur polynomials Sj{ti, . . . ,tj) are defined by the following generating relation exp(^°^-^ tjZ^) = 



Ejlo Sj{ti,t2, tj)z\ and therefore Sj = Y^l^i 



+jp=o 



■ tj^. Given a partition n = (rii, . . . , n^) G Z'j„ we 



have the Schur function Sfi{t) — <let{Sni-i+j{t))i<i.j<r- For more on the relation of these Schur functions and those 
in [2g, see [32]. 

In order to state sufficient conditions in this direction we need some preliminary definitions and results. 



Definition 3. Given a sequence C — {qI^q C 
cr e M{/S.) such that 



its Hausdorff moment problem consists in finding a measure 



Cda(C), 



i e z_i 



Moreover, if we further impose the constraint cr € A4q{A) we say that we have a restricted Hausdorff moment problem. 

Here we have made a variation in the classical definition of a Hausdorff problem, where the solutions are positive 
measures. In our Hausdorff problem we look for measures in a wider class where they do not change their sign. 
Obviously, since Mo{A) C A^(A) each solution of a restricted Hausdorff problem is also a solution of a Hausdorff 
problem. In the pages 8 and 9 in [34] J. A. Shohat and J. D. Tamarkin study Hausdorff problems and give a sufficient 
and necessary condition over the sequences to have solution. Using this result we deduce the following Lemma. 

Lemma 1. The Hausdorff moment problem for a sequence C = {ci}^Q C M has a solution if and only if 



E 

i=0 



-irc,+fc>0 V(n,fc)eZ^ 



or 



E 

i=0 



V(n,fc) e Z 



(29) 



When (291 holds a necessary and sufficient condition that ensures solution for the restricted Hausdorff moment problem 
of C is 



^ cJ < +00. 



(30) 



Proof. Theorem 1.5 in [33] states that the first set of inequalities in (29) is a necessary and sufficient condition to have 



a positive measure a solving the classical Hausdorff problem. Following their proof it is not hard to conclude that 
adding the second set of inequalities leads to a solution in A^(A). Let us take a measure a € 7W(A) and observe that 
/ i-af ^ holomorphic fimction on C \ [1, +00), then if C is its moment sequence we deduce 



Aa{t) 



= yc.x\ xe{|CI<i}. 



1=0 



/A i — Xt 

Thus, since all the c^'s have the same sign, by Lebesgue's dominated convergence Theorem we have 

da{t) 



lim 

ie[o,i) 



1 



xt 



00 



lim I x'Ci 
e[o,i) 1=0 



X 



00 

E' 

i=0 



Thus a e A^o(A) if and only if (30) takes place. 
Given the series 



□ 



00 



Wj{x) = > K,j,lX 



i=0 



X e Ai, 



,p, 



(31) 
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we introduce a set of semi-infinite matrices Qk and semi-infinite vectors Oj_k, j = • • ■ iPi k — 1, ... ,p in the following 
recursive way. First, we define 



^1,1,1 ^: 



1,1 



2,1,1 ^3,1,1 



^2,1,1 A; 



3,1,1 



-^4,1. 



^J,l 



/ -^0,^,1^ 
^1,J,1 
^2,j,l 

V : J 



j = l,...,_p. 



Then, we seek solutions 6j^2 ■= ^2,j.2 of Qi9j^2 — ^'j.i, for j — 2, ...,p, and if these solutions exist we define 



/ Ao,2,2 -^1,2,2 -^2,2,2 
Al,2,2 ^2,2,2 A3_2,2 
•^2,2,2 ^3^2, 2 A4_2,2 



Then, we look for ^-, 3 = 



A2,j,3 I which solves ^2^j,z — ^j,2i for j = 3, ... ,77, and 



3 -^1,3,3 '^2,3,3 
Al,3,3 ^2,3,3 A3, 3^3 ■■■ 
-^2,3,3 -^3^3,3 A4^3 3 



when such solutions exist we introduce 03 



In this way we get for k G {1, . . . ,p} the matrices 



V 



Qk and vectors 9j^k: j = . . . ,p, linked by Qk0j.k+i = Sj.k with expressions 



e 



k+l 



fc-|-l,fe-|-l ^l,k+l,k+l ^2,k+l,k+l ■ ■ A 
^l,k+l,k+l M,k+l,k+l ^3,k+l,k+l 
^2,k+l,k+l '^3,k+l,k+l A4_fe+i_fe+i 



V 
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Here we understand OkO 



k'^j,k+l 



\ k as 



^j,fc+i 



Ai,j,fe+i 

A2,i,fe+1 



V 



We now consider the sequences 

Ck,k '■— {Ai,A:,fe}^OJ 



M+i,k,k^i,j,k+l 

j=0 



j ^k,... ,p, 



fc = 1, . . . ,p. 



(32) 



Later, we will prove that none of the semi-infinite Hankel matrices Qk, k = l,...,p, are invertible. Hence such 
infinite linear systems are either undetermined or incompatible. In this last case we say that the systems of sequences 
{Ck,k, ■ ■ ■,Cp^k), fc = 1, . . . ,p, do not exist. 
First we need the following preliminary 



Lemma 2. The 



i=0 

converges uniformly on Aq to a function ct{x) — J diT(t)/(l — tx) corresponding to a measure a £ A^o(Ai) on Ag if 
and only if the restricted Hausdorff moment problem corresponding to the sequence {Xi : i G Z+} has a solution. 

Proof. Let us assume that the restricted Hausdorff moment problem of a sequence {A^ : i G Z_|_} has a solution. That 
means that there exists a measure a € A^o(A) such that 



A,; = / fda{t), 



i e Z+, 



lim 

xe[o,i) 



dcr(t) 

Ai 1 - 



^A, 



< +O0. 



Since |Aix'| < |Ai|, \x\ < 1 and X^i^o |Ai| < +oo, by Weirestrass' Theorem X^i^o Aj^;* converges uniformly on Aq. This 
proves the if implication in the Lemma [2] On the other hand 



lim 

.c-f 1 

a:e[0,l) 



da{t) 



l~tx 



= hm 
£ce[o,i) 



5: A. 

i=0 



because |Xi^o'^«'^N be continuous on Ag. Aj coincides with the i-th moment corresponding to the measure a 

which completes the proof. □ 
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Theorem 1. The system of weights as in (31), converges uniformly in Aq to an M-Nikishin system of 

functions if and only if for each k — 1, . . . ,p~l, there exists a system of sequences {Ck,k, • ■ • , C'p.fc) as in (32 

such that their restricted Hausdorff moment problems have solutions. 

Proof. The proof of this Theorem goes as foUows. From Lemma [2] we have that for each j = 1, . . . ,p, 

oo 

'^jA^) ^^^i.3.i^\ a; e Ao, 

1=0 

converges m Aq to a function 'sj{x) — J dsj{t)/{l — tx) if and only if the restricted Hausdorff moment problem corre- 
sponding to '■ i G 2+} has a solution. We assume that Wj^i converges uniformly on Aq to the function Sj corre- 
sponding to the Sj G A^(Ai). In order to prove the necessity in Theorem's statement we suppose that (si, . . . , Sp) = 
A4Af{ai, . . . ,ap) is an M-Nikishin system of measures as it was defined in the ^2.2 Fixed k G {1, . . . ,p} we define 
another M-Nikishin system {sk,k, ■ ■ ■ , Sk,p) = ■MAf{ak, ■ ■ ■ , cTp). Let us observe that (si^i, . . . , si_p) = (si, . . . , Sp). 

By construction for each k £ {1, . . . we have that dskj = ^k+i,jdsk,k, j = k, . . . ,p. When j = k we understand 
Sfc+i^fc = 1. Fixed j G {fc + 1, . . . Sfc+ij is a holomorphic function on C \ (— oo, 1); hence, its Taylor's series 



Wj,k+i{t) ^ h.,3,k+it\ t e Ai c Ao 



1=0 



converges uniformly to Sfc+i.j on Ai. Then, for each x £ Ai 



Sk,j[x) = yXi,j,kX ^ / Sk+i,j{t)- — = / > A,j-fc+ii— ' 



1=0 



i=0 



tx 



C30 oo „ oo oo 

— ^ ^ "^iJM+l / ^*^'dSfe^fc(t) = ^ ^ K.,j,k+lM+i,k,k 

I — n — n •'R I — n A — n 



1 = i=Q 



1=0 4=0 



which proves one implication of the equivalence. The other implication comes immediately from Lemma [2] 



□ 



We remark from the statements of Lemma [T] that the conditions in Theorem [T] are equivalent to the inequalities in 
( p9| ). Hence, by continuity criteria, such conditions are stable under perturbations of the coefficients A^^i^i, i G Z+. We 
will come to this later in §3, when we consider deformations of the weights leading to the multicomponent 2D Toda 
flows in the precise form discussed in this section. 



2.3 The Gauss— Borel factorization and multiple orthogonal polynomials 

Given a perfect combination (/i, u;i,W2) we consider [5] 

Definition 4. The Gauss-Borel factorization (also known as LU factorization) of a semi-infinite moment matrix g, 
determined by {p,,'Wi,W2), is the problem of finding the solution of 



( 



^s-^s, s = 



1 

Sifi 1 

*S'2,0 S'2,1 1 



V 
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5-1 = 



f^ofi ^0,1 ^0,2 

Sii S12 



V ; ; ; 



(33) 



In terms of these matrices we construct the polynomials 

f ^ 1. 



(34) 



where the sum is taken for a fixed a — \, 
dual polynomials 



. ,pi over those i such that a — ai{i) and i < I. We also construct the 



(35) 



where the sum is taken for a given b over those j such that b = a2{j) and j < I. 
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This factorization makes sense whenever all the principal minors of g do not vanish, i.e., if det g^'i 7^ / = 1, 2, . 
and in our case it is true because (/i,wi,W2) is a perfect combination. It can be shown that the following sets 



S 




S 



/ 'S'0,0 50,1 St) 2 
Si.i Si, 2 

S2,2 



, Si j € M, Si^ 



.0} 



are groups. Indeed, the multiplication of two arbitrary semi-infinite matrices is, in general, not well defined as it 
involves, for each coefficient of the product, a series; however if the two matrices lie on G_, the mentioned series 
collapses into a finite sum, and the same holds for G-\-. Moreover, the inverse of a matrix in 5 € G_ can be found 
to be in G_ in a recursive way: first we express 5 = 1 + X]i>o ^ii-^Y ^i^^i — diag(S'i(0), S'i(l), . . . ) a diagonal 
matrix, then we assume S^^ = I + X]i>o ^ii^Y have the same form, and finally we find that the diagonal matrix 
unknown coefficients Si are expressed in terms of ^o, . • . , 5*^ in a unique way; the same holds in G+. Given, two 
elements S S G_ and S £ G+ the coefficients of the product SS are finite sums. However, this is not the case for SS, 
where the coefficients are series. Therefore, given an LU factorizable element g = S^^S we can not ensure that g has 
an inverse, observe that in spite of the existence of 5* and S~^, the existence of S^-^S = g~^ is not ensured as this 
product involves the evaluation of series instead of finite sums. 

With the use of the coefficients of the matrices S and S we construct multiple orthogonal polynomials of mixed 
type with normalizations of type I and II 



Proposition 3. We have the following identifications 

- ^[i7i(/):£72(/-l)],a' 



A 



(0 



4(I:1l(0) 

^P2(0;^'i(i-i)],fc' 



in terms of multiple orthogonal polynomials of mixed type with two normalizations I and II, respectively. 
Proof. From the LU factorization we deduce 
I 

'^Si^igij=0, j = 0,1,. ..,1-1, Sii:=l. 

4=0 



(36) 



With the aid of (18) and (34 1 we express (36) as follows 



degA« <z.i,,(0-l, 
0<k< V2,b{^ - 1) - 1. 



(37) 



We recognize these equations as those defining a set of multiple orthogonal polynomials of mixed type as discussed 
in [2]. This fact leads to A^p^.p^^ ,^ := Aa^ where i^i = vi{l) and U2 = V2{1 — 1). Observe that for a given I each 
polynomial ^p^ji/jj.a has at much vi.ail) coefficients, and therefore we have 1171(^)1 = ^ + 1 unknowns, while we have 
|^'2(^ — 1)1 =1 equations. Moreover, from the normalization condition Sa = 1 we get that the polynomial ^[£7i;£72],ai(i) 
is monic with (iegA\^0^.j}^-^ ai(i) = f^i.ai(i)(0 ^ 1 = + 1) ~ Ij so that we are dealing with a type II normalization and 
therefore we can write 



Dual equations to ( 36 1 are 



{I) _ ^(II,ai(0) 
[!7i;i72],a 



E 

3=0 
I 

E 

3=0 



91,3 Sj^l = 1. 



1 = 0,1,. ..,1-1, 



Now, using again ( |18[ ) and ( 35 1 , ( 38 ) becomes 

[ A^^\x)vj2,hix))'Wi^aix)x''dfi{x) = 0, 



~ h=l 



deg < 1/2,6(0-1, 
0<fc<;/i,,(/-l)-l, 



(38) 
(39) 

(40) 
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while (39) reads 



(£4'^(x)«;2,b(a;))z«i,„,(z)(a;)a;'=i(')d^(x) = 1 



(41) 



b=l 



where using ( 19 ) we obtain 



fci(0 = ^i,ai(;)(' - !)■ 



(42) 




Given a definite sign finite Borel measure the corresponding set of monic orthogonal polynomials {pz}^g, degp; = 
can be viewed as a ladder of polynomials, in which to get up to a given degree one needs to ascend I steps in the 
ladder. For multiple orthogonality the situation is different as we have, instead of a chain, a multi-dimensional lattice of 
degrees. Let us consider a perfect combination (/i, wi, 'W2) and the corresponding set of multiple orthogonal polynomials 
{A[p^.p^] a}a=n '^ith degree vectors such that = |i72| + 1. There always exists compositions ^1,^2 and an integer 
I with =1 + 1 and \v2\ — I such that the polynomials {Aa-*}^!^]^ coincides with {^[i7i;i72],a}a=i- Therefore, the 
set of sets of multiple orthogonal polynomials {{^a'^''}aLii fc = 0, . . . , can be understood as a ladder leading to the 
desired set of multiple orthogonal polynomials {^[£7i;£72],o}a=i after ascending I steps in the ladder, very much in same 
style as in standard orthogonality (non multiple) setting. The ladder can be identified with the compositions {ni,n2)- 
However, by no means there is always a unique ladder to achieve this, in general there are several compositions that 
do the job. A particular ladder, which we refer to as the simplest [vi]V2\ ladder, is given by the choice ni = vi and 
n2 = 1^2 + &2,p2 ■ Many of the expressions that will be derived later on in this paper for multiple orthogonal polynomials 
and second kind functions only depend on the integers (vi, V2) and not on the particular ladder chosen, and therefore 
compositions, one uses to reach to it. 

2.4 Linear forms and multiple bi-orthogonality 

We introduce linear forms associated with multiple orthogonal polynomials as follows 

Definition 5. Strings 0/ linear forms and dual linear forms associated with multiple ortogonal polynomials and their 
duals are defined by 





(43) 



It can be immediately checked that 



Proposition 4. The linear forms and their duals, introduced in Definition^ are given by 



pi 





(44) 



a=l 



b=l 



Sometimes we use the alternative notation Q^'^ 



Q[vi:v2\ ^'^"^ Q^'"* ~ Q[v2:i'i\- It ^1^*^ trivial to check the following 



Proposition 5. The orthogonality relations 




0, < A: < z.2,&(^ - 1) - 1, 6=1,...,P2, 



0, < /c < z.i,,(; - 1) - 1, a=l,...,pi. 



(45) 



are fulfilled. 



Moreover, we have that these linear forms are bi-orthogonal 
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Proposition 6. The following multiple bi-orthogonality relations among linear forms and their duals 
hold. 

Proof. Observe that 



Q{x)Q{x) ' d/^(a;) = J S^iix)^2{x) ' S-'d^iix) 

= SgS-^ 
= 1. 



from (43) 



from (10) 



from (33) 



Definition 6. Denote by ^ , i = 1,2 the truncated vector formed with the first I components of ^i. 

We are ready to give different expressions for these hnear forms and their duals 
Proposition 7. The linear forms can be expressed in terms of the moment matrix in the following different 



1 

det 



det 



• 


1) igV+n 


)-i^[.+il 




50,0 


50,1 • • 


5o,i-i 




51,0 


51,1 •• 


5i,i-i 


f(i) 

Sl 


5i-i,o 


5i-i,i •• 


5;-i,i-i 


?1 


5i,o 


5i,i • • 


9i,i-i 





and the dual linear forms as 



/ 5o,i \ 
5i,( 

\9i-i,iJ 









VV 

/ 5o,o 5o,i 
51,0 51,1 



1 



det 



det 



5o,i-i 
5i,i-i 



5i-i,o 5i-i,i ■ ■ ■ 5;-i,i-i 



fC-i) 

?2 



5o,i \ 
5i,i 



5/-1,/ 



I > 1, 



/ > 0. 



Proof. See Appendix \K\ 

As a consequence we get different expressions for the muhiple orthogonal polynomials and their duals 
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Corollary 1. The multiple orthogonal polynomials and their duals have the following alternative expressions 



and 



A 



(0 





[gift 


giA ■ ■ ■ giA-i 




Xl,a 


= Su (0 





.. 


-1 ['+1] 








On n On i • • • 

yu,u i/u,l 


On 1 1 

yu,(— 1 


(0) 
Al.a 


1 

~ dct 




51,0 51,1 ■■• 






det 


yt— 1,0 y/— 1,1 


yt— 1 ,i— 1 


(i-i) 






\ gifi gi,i ■ ■ ■ 


gid^i 
( go.i \ 


(0 

Xl.ai 






- ix^Vig^'^r' 


51,/ 


) 



\gi^i,i) 







det 



det 




VV 

/ 5o,o 5o,i ■ ■ • 5o,i-i 

51,0 51,1 ••• 5i,/-i 

5i-i,o 5i-i,i • • • gi-id-i 



JO) (1) 
\ A2,f) X2M 



go,i \ 
gi,i 



5/-1,/ 



Observe that (208), Appendix [A] implies 



5, 



i.i 



(/-I) 

X2,6 



detg['+il 
detgW 



(0 

^2,6 y 



2.5 Functions of the second kind 



/ > 1. 



(49) 



(50) 



I > 0. 



(51) 



(52) 



The Cauchy transforms of the linear forms ( 44 ) play a crucial role in the Riemann-Hilbert problem associated with 
the multiple orthogonal polynomials of mixed type [14j . Following the approach of Adler and van Moerbeke we will 
show that these Cauchy transforms are also related to the LU factorization considered in this paper. 

Observe that the construction of multiple orthogonal polynomials performed so far is synthesized in the following 
strings of multiple orthogonal polynomials and their duals 



V ■■ J 



■ = 



(1) 



V J 



= {S'^)^X2,b, a^l,...,pi, &=l,...,p2- 



(53) 



In order to complete these formulae and in terms of x* as in ([8]) we consider 
Definition 7. Let us introduce the following formal semi-infinite vectors 



'£h — 



c, 



(1) 



V ■■ J 

that we call strings of second kind functions. 



(1) 



V ■■ J 



= [S ^VxiAz), b=l,...,p2, a = l,...,pi, 



(54) 
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These objects are actually Cauchy transforms of the linear forms Q^^\ I G Z+, whenever the series converge and 
outside the support of the measures involved. Notice that fixed z E C the entries in each string '^a and "^^b are series 
not necessarily convergent. In the non-convergent case we obviously understand the definition only formally. For 
each I € Z+ we denote by Z)a ' and D^'"* on C the domains where the series C« and Ci'^ are uniform convergent, 
respectively, and we understand them as their corresponding limits. From properties of Taylor's series, we know that 
uniform convergence of these series hops only on _Da ■* and Z?^'^ when they are the biggest open disks around z = oo 
which do not contain the respectively supports, supp(w2,ad/x) and supp(u>2,6d/i). Outside the sets Da'^ and D^l^ the 
series diverges at every point. Hence to have non-empty sets in _Da ■* and ZJ^'' the corresponding supports supp(i(J2,ad/J.) 
and supp('u;2,fed/z) must be bounded. 

Proposition 8. For each I G 2+ the second kind functions can be expressed as follows 



Z — X 



Proof. The Gauss-Borel factorization leads to 

oo / 



d/i(a;), z e d'^P \ supp(wi,6d^(x)), 
d^(a;), z e D^^ \ swpp{w2,ad^i{x)). 



(55) 



n=0 fc=0 



= E / E'^'fe^'''^'''^i.«i('=)(^)"'2,fc(2;)^;;TTdM(a;) 

71=0 k=0 

E^ / X^^Q^'\X)W2A^)M^)- 

n ^ 



use ( 12 1 



use (441 



n=0 



When D^f^ \supp(w2,6d/i) = the proof is trivial. Given a non empty compact set /C C \supp(w2.hd/i) ^ and re- 
calling the closed character of supp(w2.hd/i), we have that the distance between them d^'"' {K) :— distance(/C, supp(w2,6d/i)) > 
is positive and that sup{|z| : z G /C} =: M/c < +oo. Taking into account that the series 

oo 

. n ^ J 



(0^ 



converges uniformly on K, we can ensure 



lim sup 

|z|G/C 



x''Q^'\x)w2,b{x)d^i{x) 



= 0. 



(56) 



Hence, we have the bound 



E^ J x'Q^'Hx)w2A^)Mx) ~ I Q'^'Hx)w2,b{^)^^Mx) 

x''Q^'Hx)w2,bix)d^iix) 



/ .»gW(,)^,,,(,)M?i) 



< 



< —77-. sup 



4\JC) \.\eK 



,\fz G /C. 



(57) 



Taking into account (56) we deduce from (57) the first equality for any compact set K.. Therefore, we get the first 
claim of the Proposition; the second equality can be proved analogously. □ 

Given I > 1 and a — 1, • • • ,p the + (— ) associated integer is the smallest (largest) integer l^a (l-a) such that 
> I {l-a < I ) and a(Z+a) = a (a(/_a) — a). It can be shown that 



a < a{l), 

l-a ■= {l, a = a(/), 

,9(01^1 -ELa-t-i^i - 1' a>a(/-l), 
'm + l)\n\+J:^Zln^, a<a{l), 
l+a ■= {I, a = a{l), 



(58) 
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To give a determinantal expression for these second kind formal series we need 
Definition 8. We introduce 



^k,a— 9k' MZ Oa,{k'),a, ^ k,b — 9k,k' Z d, 



a-2{k'),b- 



(59) 



k'=l. 



Here l-^-a is the + associated integer within the ni composition, while l+b is the + associated integer for the n2 
composition. 

With these definitions we can state 

Proposition 9. The following determinantal expressions for the functions of the second kind hold 





^ 90.0 


90,1 ■ 


■ 9o,i-i 


f(0 

0,6 


rpr det 

det 


91,0 


91,1 ■ 


■ 


1,6 


9i~i,o 


9i-i,i ■ 


• 9i-i,i-i 


^ Z-1,6 




y 9lo 


9i,i ■ 


■ 9i,i~i 


f(0 

^ l,b' 



I > 1, 



(60) 





/ ffo.o go,i ■ 


• 5o,/-i 


90, 1 


rrm det 

det 


91,0 gi,i ■ 


• 9i,i-i 


gi,i 


gi-1,0 gi-1,1 ■ 


■ 91^1,1-1 


9i-i,i 




y ^ 0,a ^ l,a ■ 


' ^ l-l,a 


^ l,a 



I > 1. 



Proof. See Appendix [A] 

Following [2^ we consider the Markov-Stieltjes functions and polynomials of the second type. 
Definition 9. The Markov-Stieltjes functions are defined by 



(61) 



□ 



fl-a.b{z) 



Z — X 



-d^{x), 



(62) 



in terms of which we define 



HI'\z) ■.= Y.^aHz)f^a.b{z)-Ci'\ 



a=l 

P2 



(63) 



6=1 



Proposition 10. The functions h9'^ and Ha'^ are polynomials 



Proof. The reader should notice that the functions fff'"' and Ha' are 



Hi'\z)= /£™l,a(x) 
a=l 



A^I>iz)-A'^>ix) 



(I), 



-W2,b{x)dfi{x), H'^a\z) = / Wi,a(x) 

6=1 



A'^i\z)~A^^>{x) 



(0/ 



W2,b{x)<:l^{x), 



and as = a; is a zero of the polynomials Aa (z) — Aa {x) and A\^ (z) — A\^ {x) from the above formulae we conclude 
that they are indeed polynomials in z. □ 
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2.6 Recursion relations 

The moment matrix has a Hankel type symmetry that imphes the recursion relations and the Christoffel-Darboux 
formula. We consider the shift operators defined by 

oo 
fe=0 

notice that 

• A(j leaves invariant the subspaces Ila/IR"", for a' ~ I, . . . ,p, and Ha'^a = Aalla'. 

• The set of semi-infinite matrices {A^}a=i p is commutative. 

i=i,2,.'.. 

• We have the eigenvalue property 

AqXo' = 5a,a'ZXa- (65) 

Definition 10. We define the following multiple shift matrices 

Pl P2 

Ti:=^Ai,a, T2:=^A2,b, (66) 

a=l 6=1 

and we also introduce the integers 

Ni^a ■= \ni\ - ni^a + I = ?^i.a' + l, a = l,...,pi, iVi := nrax Ni^a, 

a' = l,...,pi 

a' 

N2,b ■= \n2\ - n2.b + 1 = »T-2.fc' + 1, & = 1,---,P2, := max N2,b- 

b =l,...,p2 
b'^b 

A careful but straightforward computation leads to 
Proposition 11. We have the following structure for Ti and T2 

Ti = Di^oA + i?i,iA^i ^ + • • • + i?i,p, A^i-i , T2 = D2,oA + D2sA^'-' +■■■ + D2,p,K^''^^ . 

where -Di,a, a = 1, . . . and D2.b, b — \^ . . . ,p2, are the following semi-infinite diagonal matrices: 

7^i,=diag(i?uo),^i,a(i),...), DUn)---\i' ":^!l^!tS'='"'"'~!' f^^' ^i.o=i-i:/?i,a, 

£'2,b = diag(i:»2,&(0),I?2,b(l),---), D2,b{n):^< I ' £'2,0 = I - > -D2,fc- 

1^0, n ^ /c|n2| + }J6' = l"2,6' - 1, fceZ+, ^ 

In terms of these shift matrices we can describe the particular Hankel symmetries for the moment matrix 
Proposition 12. The moment matrix g satisfies the Hankel type symmetry 

T,g = gTj. (67) 



Proof With the use of ([19|) and ([12|) we get 

Ai,agn2,b = Ui.agAl^, (68) 

and summing up in a = 1, . . . ,pi and b — 1, . . . ,p2 we get the desired result. □ 

Observe that from (67) we deduce that in spite of being all the truncated moment matrices gW, I — 1,2,... 
invertible, the moment matrix g — lim/_>.oo ff^'^is not invertible. Suppose that the inverse g~^ = (gi,j)i,j=o,i,... of g 
exists so that (67) implies g~^Ti — Tj g~^, and therefore gi^ — goj = for all i, j — 0,1, . . . , which is contradictory 
with the invertibility of g. 
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Proposition 13. From the symmetry of the moment matrix one derives 

S'TiS'-i = STjS-'^. (69) 



Proof. If we introduce ( 33 ) into ( 67 ) we get 

□ 

Definition 11. We define the matrices 

J:=J+ + J-, J+ := {STiS-^)+, J_ ;= (5tJS*-1)_, 

where the sub-indices + and — denotes the upper triangular and strictly lower triangular projections. 

Thus, J-f is an upper triangular matrix and a strictly lower triangular matrix. Moreover, from the string 
equation ( 69 ) we have the alternative expressions 

We now analyze the structure of J+ := {SDi^o^S-^)+ + {SDi^iA^^-^ S-'^)+ + ■■■ + {SDup.A^^^pi S-'^)+. It is clear 
that we need to evaluate expressions of the form SEijS~^ with i = Ki(k,a) — 1 and j ~ Ki(k + 1, a — 1) being 
Ki(fc, a) := k\ni\ + J21'=i '^i,a' ■ Given the form of S, see (33), we have 

for some numbers si^i' € M depending on the coefficients of S and on i, j; this matrix has zeroes everywhere but on a 
region of it that can be represented as a right triangle with hypotenuse lying on the main diagonal, this hypotenuse 
has its opposite vertex precisely on the {i,j) position. Therefore 

Pi oo / \ 
a=l k=0 \ ld'ehi,k,a / 

We see that J+ can be schematically represented as a staircase, the ni-staircase, descending over the main diagonal 

with steps -which are built with right triangles with hypotenuse lying on the main diagonal and opposite vertex 

(and therefore corner of the step) located at the (Ki(fc, a) — 1, Ki(fc + 1, a — 1)) position of the matrix- having width 

and height given by the integers in the composition fii. For example, the j-th step has width j__rj_i and height 

' PI ^ PI J 

j+i_rj+ii. A similar description holds for jj but replacing the composition ni by n2. Therefore, the matrix J 

'PI PI ' 

is a generalized Jacobi matrix and, in contrast with the non multiple case, now is multi-diagonal (having in general 
more than three diagonals) and has a diagonal band of length A^i + A'2 + 1. Moreover, this band has a number of 
zeroes on it, according to the ni-stair on the upper part and to the ri2-stair in the lower part, we refer to this as a 
double (ni, n2)-staircase shape. To illustrate this snake shape let us write for the case Hi = (4,3,2) and 7x2 = (3,2) 
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the corresponding truncated, I — 27, Jacobi type matrix 



Jl27] 
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where * denotes a non-necessarily null real number. We can write 

J = Jjv, A^i + • • • + JiA + Jo + J_iA^ + • • • + J-w, (A ' ) 



where Ji — diag(Ji(0), . . . ). For convenience we extend the notation with Jr{s) — whenever r + s < or 

We introduce 

Definition 12. The semi-infinite vectors Cb and Ca are given by 

Cb ■■= Se2.b{0), 6=l,...,p2, 
Ca (S'~^)^ei,a(0), a = l,...,pi. 

It is not difficult to show that 

ri2,lH l-n2.b — 1 "l.lH hiii.a — 1 



Cb 



<5';,n2_i+---+ri2,i,-iei, Ca — ^ 

;=o 1=0 
The semi-infinite matrices J and have the following important property 
Proposition 14. The following equations are fulfilled 

j£/a{z) = Z£/a(z), S>/b{z) = ZS^b{z), 

J'^biz) = Z<^b(z) - Cb, J'^'ia^z) = Z'^aiz) " 



(5-1) 



"l.lH hni, 



-1,/e/. 



Proof. From (53) and (54) 



J^aiz) = STiS ^SxiA^) zSxi^a = Z£/aiz), 

J%{z) = STjS-'Sx*2,b{z) = S{zxlb{z) - ebiO)) = z%{z) - Cb, 
where we have taken into account that Tjx2 bi"^) ^ -^^2 bi^) ~ Cb{0). For we proceed similarly: 

J^M,{Z) = {S-Y^2S^{S-Yx2.b{z) = z{S~Yx2,a - Z^b(z), 

J'^^aiz) = (S-yrJS'^iS-Yxlaiz) - iS-Yizxla{z) - CaW) = Z<#a(z) - C^. 
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Theorem 2. The multiple orthogonal polynomials and their associated second kind functions fulfill the following 
recursion relations 



(/-JV2) 



(^) 



while the dual relations are 



+ J^.,(/M^^^)(z), 



z&J\z) -ca^ J-nAI + N2)C'i+^^-\z) + . . . + Jjv,(/ - N,)Cf'\z)- 



(75) 



(76) 



We see that given integers {Vi, V2) there are several recursion relations associated with A^p^^.p^j^a- In fact they are as 
many as different ladders exists leading to this set of degrees. For the simplest ladder, i.e. fii = vi and 7^2 = ^'2 + e2,p2; 
we get the longest recursion, in the sense that we have more polynomials contributing in the recursion relation, 
as smaller are the integers in the compositions shorter is the recursion. Observe also that the multiple orthogonal 
polynomials involved in each case are different. 



Attending to (70) we get that the recursion relations corresponding to ^ = 8 an / = 14 are of the form 

a = 1,2,3, 
a = 1,2,3. 



zA^^\z) = *AW(z) + . . . + *Aii5)(z) + A(i6)(,), 
zA^^{z)^.A^{z) + --- + .A^^{z), 



We see that the first recursion has 13 terms while the second one only 7 terms. 

In order to identify these polynomials with mops of the form A^^^^.p^^^] we use following the table of degrees for the 
compositions Hi = (4, 3, 2) and n2 = (3, 2) is 



1 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


02(1 - 1) 


(4,0,0) 
(2,0) 


(4,1,0) 
(3.0) 


(4,2.0) 
(3,1) 


(4,3,0) 
(3,2) 


(4.3,1) 
(4,2) 


(4,3,2) 
(5.2) 


(5,3,2) 
(6.2) 


(6,3.2) 
(6,3) 


(7,3.2) 
(6.4) 


(8,3,2) 
(7,4) 


(8,4,2) 
(8,4) 


(8.5,2) 
(9.4) 


(8,6.2) 
(9,5) 


(8,6.3) 
(9.6) 


(8.6,4) 
(10,6) 



2.7 Christoffel— Darboux type formulae 



From ( 53 1 and ( 54 ) wc can infer the value of the following series constructed in terms of multiple orthogonal polynomials 
and corresponding functions of the second kind. 

Proposition 15. The following relations hold 

00 

5:cW(z)4)(z') = f^, \z'\<\z\, 

^Cl"(z)A«(.') = ^, \z'\<\zl (77) 



1=0 



1=0 



^C(')(z)C«(z') = 1^1^ |,,| > ^^^^^ 



where Ra.b is the radius of any origin centered disk containing supp('U(i^a'U(ci^;,d/x). 

Proof. See Appendix \K\ □ 

2.7.1 Projection operators and the Christoffel— Darboux kernel 

To introduce the Christoffel-Darboux kernel we need 
Definition 13. We will use the following spans 

= Kur\-..,e^'^}, ^mf\---,it'^}. (78) 

and their limits 

^i = { E Q^f'.Qeia}, ^2 = { 5] Q#:C,eM}. (79) 

0<i<Soo 0<;<goo 
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The corresponding splittings 

ni=Hf^ ®{HfY, H2=Hfl®(HfV, (80) 

induce the associated orthogonal projections 

tt['^ : Hi ^ H^^ , (81) 

In the previous definition I <ti oo means that in the series there are only a finite number of nonzero contributions. 
It is easy to realize that 

Proposition 16. We have the following characterization of the previous linear subspaces 

h'') =M{Q(°),...,g('-i)}, =M{Q("',...,Q('-i)}, 



o<;<oo o<;<goo 
Definition 14. The Christoffel-Darboux kernel is 

i-i 



fc=o 

This is the kernel of the integral representation of the projections introduced in Definition [TSj 
Proposition 17. The integral representation 

i4^f){y) = / K^Hx,y)f{x)dfi{x), V/ e Hi, 



{4^f)iy) = / ifW(2/,x)/(x)dA.(x), V/ e n2 



holds. 



(82) 



Hi = { E c/Q('\QeM}, ^2 = { E c,Q('\QeM}. (83) 



ifW(x,2/):=E^?^'^(y)<3^'^(^)- (84) 



(85) 



Proof. It follows from the bi-orthogonality condition ( 46 1 . □ 

This Christoffel-Darboux kernel has the reproducing property 
Proposition 18. The kernel A'W (x, y) fulfills 

K^^Kx,y)= J K^^\x,v)K^^^{v,y)dn{v). (86) 

Proof. From 

f{y) = I K^Hx:y).f{x)dfi{x), V/ e nf\ 

f(y) = j K^'\y,x)f{x)df,{x), V/ G H^^ , 

&ndK^^^{x,y) e Tif as a function of y and if ^ (cc, y) € T-Ly as a function of x wc conclude the reproducing property. □ 
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2.7.2 The ABC type theorem 

We also have an ABC (Aitken-Berg-CoUar) type theorem — here we foUow ^35l — for the Christoffel-Darboux kernel 
Definition 15. The partial Christoffel-Darboux kernels are defined by 

i-i 



(87) 



fc=0 



Observe that 



a=l,...,pi 
b=l,...,P2 

We introduce the notation 
Definition 16. Any semi-infinite vector v can be written in block form as follows 

( „W \ 



,[>'] 



(89) 



where ut'l is the finite vector formed with the first I coefficients of v and ut-'l the semi-infinite vector formed with the 
remaining coefficients. This decomposition induces the following block structure for any semi-infinite matrix. 



(90) 





g[l>l] 


gl>l.l] 


g[>l] 



Prom ( 33 ) we get 

Proposition 19. Given a moment matrix g satisfying ( |33[ ) we have 

and (5-i)W = (5W)-\ (5-i)[^'l = (5[^'l)-i. 
Proof. Use the block structure of g, S and S. 

Then, we are able to conclude the following result 
Theorem 3. The Christoffel-Darboux kernel is related to the moment matrix in the following way 

Proof. The ABC theorem is a consequence of the following chain of identities 



(91) 



□ 



(92) 



<L(a;,y) = (n['W,(:.))T(nW<(y)) 

= xld^) (n w 5- In w ) (n w sui'^)xiAy) 
= ixl\{AV{9^'Y'x^^Jy) 



the sum is over the first I components 



see (|43l 

lower and upper form of S and S 



LU factorization ( 33 1 . 



□ 



We immediately deduce the 
Corollary 2. For the Christoffel-Darboux kernel we have 

iff"(^,2/) = (4'k^))^(3"i)~'el"(y)- 



(93) 
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2.7.3 Christoffel Darboux formula 



In this subsection we derive a Christoffel-Darboux type formula from the symmetry property (67) of the moment 
matrix g. We need some prehminary lemmas 

Lemma 3. The relations 

_ (TW)T(gW)-i = (g['l)-i(5['->'l(TP")^ - tP'15[^'''1)(5W)-\ (94) 

hold true. 



Proof. The first block of (67) is 

xWgW + t['->'1^[>/,;] ^ ^m(TW)T + ^[i,>i](y[;,>ii)T^ 

from where the result follows immediately. 
Lemma 4. We have 

Proof. It follows from the block decomposition of Definitions [16] and the eigen- value property of Tg 



□ 

(95) 
□ 



After a careful computation from Definition 10 we get 
Lemma 5. If we assume that I > max(|ni|, \n2\) we can write 

pi 



^2 - Z^^(/-l)_,%6 



(96) 



b=l 



Here l±a is the ± associated integer within the ni composition, while l±b is the ± associated integer for the n2 
composition. 

Finally, to derive a Christoffel-Darboux type formula we need the following objects 
Definition 17. Associated polynomials are given by 



(<?['i)-^x^(y) 



(x^'.;/'(a^))'(.9"+^0-^e,_„, 



1(0 



( 9o.i+, \ 

\.9/-i,r, 6/ 



(97) 



with the corresponding linear forms given by 



P2 



a' = l 
Pi 



b' = l 
P2 



(98) 



a' = l 6' = 1 

Then, we can show that 
Theorem 4. Whenever I > max(|ni|, \n2\) the following Christoffel-Darboux type formulae 

P2 Pi 



(x - y)i^f 1 (x, y) = E ^+M' (^)^-M' - E ^- J (^)4la' iy) > 

6=1 a=l 

P2 Pi 

(x - y) = E (y) - E '^(^)Q+i(2/)' 



6=1 



o=l 



hold. 
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Proof. From Lemma [3] we deduce 

(xf;.^(-))^((5W)-^T['l-(Tfb^(.gW)-i)xt(y) = (X^^ 
so that, recalling Theorem [Sj we get 



(99) 



or 



Finally, from Lemma [5] we conclude 

pi 

a=l 

Pi 

a=l 

P2 



5/ + a,/-l) 



6=1 



5 



b=l 



9iJ+b 



and consequently 



b=l 



( 9o,h, \ 
9iJ+b 

\9i-i.i.i/ 



pi 



a=l 

(100) 



Recalling Definition [TT] we get the announced result. 

The associated linear forms are identified with linear forms of multiple orthogonal polynomials as follows 
Proposition 20. We have the formulae 



□ 



^(0 _ ^(II,a) ^(0 _ ^(I,h) ^(0 _ ^(II,&) Ml) _ n^ha) 



!+a — ^[,7i(;_i)+g>i^;£72(i-l)]' ^-f, — '^[i7i(0;i72(0-e2,b] ' ^ +b ~ p2 (i- l)+e2,6 (I- 1)] ' 



(101) 

Proof. See Appendix [A] □ 
Proposition |20] allows us to give the following form of the ChristofFel-Darboux formula stated in Theorem |4] 
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Proposition 21. For I > niax(|rii|, \n2\) the following 

P2 



{x y)K^^^{x,y) ~^Qip^(Li)+e2,biPi(i-i)]^^^QlpJi-i)iff2{i-i)-S2,b]^y^ 



6=1 



X!*5[£72(i-l);£7i(i-l)-ei,„](2^)'3[^'i(i-l)+ei.„;i72(i-l)]'^^)- 
a=l 

(x - y)K\)y{x,y) = ^\a^(i-i-)+e^ ^-.v, a-ii].b' (^)^p', »-ii:i7,»-iWeo^].a' (y) 
Pi 

/ . "^[i72(i-l):fl(/-l)-ei „] ,h' '^"^''"^pi (i- l)+ei c\V2{l-l)],a'^y > ' 



P2 



f)=l 



(102) 



(103) 



holds. 



Relation (102) is precisely the Christoffel-Darboux formula derived in [T3], the difference here is that derivation 
is based on the Gauss-Borel factorization problem for the moment matrix; i.e. only on algebraic arguments, and 
not in the Riemann-Hilbert problem found in [IT, and hence the conditions on the weights are not so restrictive. 
However, the reader should notice that the Christoffel-Darboux kernel does not depend on the ladder determined by 
the composition vectors ni,n2, but only on the degree vectors vi{l — 1) and ^2(1 ~ !)• This was noticed in |15| for 
type I multiple orthogonality. 

Proposition 22. The associated polynomials introduced in Definition \l'l/\ have the following determinantal expressions 



( 



A 



(0 



det gV\ 



det 



ffo.o ffo.i • • • ffo,;-i 
ffi.o 51,1 ■•• 

9i-ifi 91-1,1 ■ ■ ■ 9i-id-i 



\ 9l + afi 9l + a,l • • • 9l + a,l-l 



(0) \ 

1 r,' 



XI. a 



(i-i) 



(104) 



(-1) 



det gl'+i] 



det 



50,0 ■ 


■ 90.1-1 


90.1 


9i-a~i.o ■ 


■ 9l-a-l,l-l 


9l-a-l,l 


ffi-a + 1,0 • 


■ 9l-a + lJ-l 


9l-a + ld 


ff;-i,o • 


■ 91-lJ-l 


91-1,1 


(0) 
X2,b' ■ 


(l-l) 

■ X2M' 


(0 



(105) 





' 90.0 ■ 


1 

■ 9o,i-b-i 1 5o,r_6+i ' 


■ 90,1-1 


xt' \ 


" det^I'+U '^^^ 


91-1,0 ■ 


1 

: 1 : 

i 

' 9i-i,T_t-i \9i-i,i_h+i ■ 


■ 91-1,1-1 


xtP 




\ 91,0 ■ 


1 

■ 9i,i_t-i 1 5iJ_i,+i ■ 


■ 91,1-1 


xS. I 



(106) 



/ 5o,o 
ffi.o 



A 



(0 1 
+'''^' detgW 



det 



90,1 
91,1 



90,1-1 
91,1-1 



91-1,0 91-1,1 ■ ■ ■ 91-1,1-1 



AO) 



(1) 
A2,6' 



(l-l) 
X2.V 



9o,i+b \ 
9i,i+b 

9i-i,L, 



X2.b' 1 



(107) 
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3 Connection with the multi-component 2D Toda Lattice hierarchy 



In this section we introduce deformations of the Gauss-Borel factorization problem that give the connection with the 
theory of the integrable hierarchies of 2D Toda lattice type, in the multi-component flavor case. First, we introduce 
the continuous flows and then the discrete ones. Let us stress that both flows could be considered simultaneously but 
we consider them separately for the sake of simplicity and clearness in the exposition. 

3.1 Continuous deformations of the moment matrix 
Definition 18. The deformed moment matrix is given by 

9n,,nAt) ■■= Wofi,{t)gWo,nM"\ (108) 
where we use the following semi-infinite matrices 

Pi OO P2 CXD 

a=l 3 = 1 b=l 3 = 1 

depending on t — {tj^a, ij^b)j.a.b with tj,a, ij.b G li; J = 1, 2, . . . , a = 1, . . . ,pi and 6=1,... ,p2. 

As in the previous section and when the context is clear enough we will drop the subscripts associated with the 
compositions fii and n2. The reader should notice that the following semi-infinite matrices are well defined 

Pl OO P2 OO 

a=l 3 = 1 b=l 3 = 1 

This deformation preserves the structure that characterizes g as a moment matrix, in fact we have 
Theorem 5. The matrix g{t) is a moment matrix with new "deformed weights" given by 



OO 

Wi^a{x,t) = Sa{x,t)wi^a{x), Sa := CXp ( ^j^gX^ ) , 

OO 

W2,b{x,t) = <ib{x,t)~'^W2,b{x), ib := exp ['^tj^bX^). 



(109) 



Proof. Observe that 

Pl P2 

j">Oa=l j>Ob=l 

where trj"' is the j-th elementary Schur polynomial in the variables tja and a\^^ is also an elementary Schur polynomial 
but now in the variables —tj i,. To prove (109) we first discuss the action of Ai ^ and A^j^ on g explicitly. Recalling 
(12) it is straightforward to see that 

(Ai,agA]^b)ij = J a;'''^'^+^Wi,ai(»)(a;)w;2,a2(j)(a;)a;''"'^'^+Ma,(,)^a'5a2 0).()dAi(a;), 

and consequently the following expression holds 

Pl P2 „ 

iWo9W,-\, = EE / ^'^'^^(E^^"^^0^i''^^«('^)^2''^^(^)('^K E ^^I^x"^)x''^^^^S,^^,)^J,^^^)^,dKx), 

a=l 6=1 i>0 m>0 



that leads directly to (109). 

□ 
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That the sign definition of the weights is preserved under deformations is ensured by the fact that all times t 
are real. Let us comment that these deformations could be also be considered as evolutions, and from hereon we 
indistinctly talk about deformation/evolution. If the initial measures have bounded support then there is no problem 
with the exponential behavior at oo of the <ff factors; however, for unbounded situations a discussion is needed for 
each case. 

The Gauss-Borel factorization problem 

gn,,nAt)^Sit)-^S{t), G G_, S{t) e G+, (110) 

with S(t) lower triangular and S{t) upper triangular, will give the connection with integrable systems of Toda type. 



Let us assume that the weights in (wi,W2) are of the form (28) and that conform an M-Nikishin, then Theo- 
rem [T] indicates that for small values of the times the new weights are also in the M-Nikishin class, ensuring that 
{'Wi{t) , W2{t) , fi) is a perfect system and therefore the Gauss-Borel factorization makes sense. 

3.2 Lax equations and the integrable hierarchy 

Let us introduce the Lax machinery associated with the Gauss-Borel factorization that will lead to a mult i- component 
2D Toda lattice hierarchy as described in [28 : 

Definition 19. Associated with the deformed Gauss-Borel factorization we consider 

1. Wave semi-infinite matrices 

W{t) S{t)Wo{t), W{t) S{t)Wo{t). (Ill) 

2. Wave 

*,(z, t) W{t)xiAz). ^b{z, t) := Wit)xlbiz), (112) 

and adjoint wave semi-infinite vector function^ 

^:{z,t):^iW{t)-Yxlaiz), niz,t):^{W{t)-Yx2A^)- (US) 



3. Lax semi-infinite matrices 

La(t) := S{t)ki^aS{t)-\ hit) := S{t)Al,S{t)-\ (114) 

4-. Zakharov-Shabat semi-infinite matrices 

Bj,a := {Li)+, B,., (L^,)_, (115) 

where the subindex -\- indicates the projection in the upper triangular matrices while the subindex — the projection 
in the strictly lower triangular matrices. 

Observe that 

La^a' = Sa,a'Z<fa', LJ = St^b^Z-'^l,. (116) 

We also mention that the matrices S and S correspond to the Sato operators (also known as gauge operators) of the 
integrable hierarchy we are deling with. Some times |^ the operators La are referred as resolvents and the Lax name 
is reserved only for a convenient linear combination of the resolvents. 

The reader should notice that as S{t) G G_ and Wo{t) e G+ the product W{t) = S{t)Wo{t) is well defined as 
its coefficients are finite sums instead of series, for {W{t)~^)^ = {S{t)~^)^ {Wo{t))~^)^ we can apply the previous 
argument and therefore the product is well defined. However, {W{t)~^)^ = {S{t)~-^)^ {WQ{t)~-^)^ is a product of 
elements which involves series instead of finite sum and its existence is not in principle ensured. The situation is 

^In this point the reader should notice that there are two differences between this definition of wave functions (also known as Baker- 
Akheizer functions) and the one common in the hterature, see for example |4]. Our modifications are motivated by two facts, i) we prefer 
to be a polynomial in z and not in z~^, up to plane-wave factors, ii) we choose to have a direct connection between wave functions 
and Cauchy transforms of polynomials, with no factors multiplying the Cauchy transforms when identified with wave functions. If we 
denote by small i/i the wave functions corresponding to the scheme of for example [1] then we should have the following correspondence 
« i'i'Hz), z{^l)m{z) « « i'i'Hz) and Cl'*)(0(2-l) « (^*)(0(^). 
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reproduced with W{t) — S{t)Wo{t), and the existence of the product is not guaranteed. However, we notice that 
the simultaneous consideration of the factorization problems (110) and (33) leads S{t)^^S{t) = Wo{t)S^^ SWo{t)^^ 
that shows two products involving series, namely Wo{t)S~^ and SWo{t)^\hut they are well defined if we assume the 
existence of both LU factorizations. From hereon we give for granted the existence of W and W~^, and as we will see 
they indeed involve series, which in the convergent situation lead to Cauchy transforms. 



Proposition 23. For the wave functions we have 
vI/«(z,t) = AW(z,iK(z,t), 



(117) 



where A^a'^x , t) , A''jI'\x , t) are the multiple orthogonal polynomials and dual polynomials (in the x variable) correspond- 
ing to (109). The evolved linear forms, associated with weights (109), are 



pi 



a— 1 a— 1 

P2 P2 

Q^^Hx^t) Y.^Al)^'^{x,t)w2,b{^,t) = E(*^)^'^(^'*)^2,6(x) 

b=l 

which are bi-orthogonal polynomials of mixed type for each t 

Q^'\t,x)Q^''\t,x)Ati{x)^5i^k, 



b=l 



l,k>0, 



and 



¥,'\z,t) 



QW(a;,t) 



W2,bix)dn{x), 



(*:)W(z,t) 



Wi^aix)dn{x). 



Proof. From the definitions ([112) and (113), and the factorization problem Wg = we conclude 



*6 = Wxlb = SiWog)xh^ 
We get, in terms of the linear forms, the following identities 



(118) 
(119) 

(120) 

(121) 
(122) 



¥,'\z,t) = 



QW(x,t) 



W2,b{x)dfJ.{x), 



(*:)«(z,t) = 



QW(x,i) 



Wi^a{x)dfl{x), 



where the Cauchy transforms are understood as beforej^ 



□ 



We must stress in this point that these functions are not the evolved second kind functions of the linear forms 



Ci'\z,t) 



Q<^^\x,t) 



Z — X 



W2,hix,t)dfl{x) 



gW(x,t) 



wi,a{x,t)d^{x). (123) 



Theorem 6. For = 1, 2, . . . , a, a' 1, . . . ,pi and b,b' — 1, . . . ,p2 the following differential relations hold 



^The reader should notice that there is a difTcrcnce in this semi-infinite context, appropriate for the construction of multiple orthogonal 
polynomials, and the bi-infinite case which is the one considered in 1371 . In the present context wc do not have expressions, as we do have 
in the bi-infinite situation, of the form 



%''>{z,t) = (Po-f Piz-i -f ■••)exp(^ij-6^^), 

{z, t) = {Qo + Qiz-^ + ---)e^p{~J2 • 
The reason for this issue is rooted into non-invertibility of A^. Indeed, for the semi- infinite case, we have 



exp ( Cj{Ajy)x*a = [cxp ( CjZ^)xl] 

j=l j=l 



where the subindex — stands for the negative powers in z in the Laurent expansion; while in the bi-infinite case we drop the — subindex 
in the previous formulae. 
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1. Auxiliary linear systems for the wave matrices 



dW 

dtj.a 



dW 



dt 



dW 



Bj.aW, 



2. Linear systems for the wave and adjoint wave semi-infinite matrices 



dtj^a 
9t:,,a 



-^7, a ^ a 5 



J, a a'' 



dtj^b 
dtj^b 



j^b ^ a' : 



],b ^ a' ' 



d'^b' 

dtj.a 

dtj„ 



= B.a'^b', 



dW 
9tj,b 



dtj^b 



dt 



Bj,bW. 



-Bibn- 



3. Linear systems for multiple orthogonal polynomials and their duals 



4-. Lax equations 

dLa' 



{Bj,b)K' 



dt 



-B 



d£/b' 



0,0. 



dt 



3,b 



dtj^a 



— [Bj,a-, La'], 



dtj,b 



~ [Bj,b,La'], 



dLb- 

dtj.a 



— [BjM, Lb'], 



dLb' 



dt 



]-b 



— [Bj,b, Lb']- 



5. Zakharov-Shabat equations 



dB,,a 


dB,',a' 


dtj'^a' 


dt],a 


dBj^b 


dBj'.b' 


dify 


dtj,b 


dBj,a 


dSj'^b' 


dtj'^b' 


dtj,a 



[Bj,a, Bj'^a'] — 0, 

- [Bj.b, Bj'^b'] — 0, 

[Bj,a,Bj'^b'] = 0. 



Proof. To prove ( 124 1 we proceed as follows. In the first place we compute 

dWo 



dti,a 



dt 



O-b 



i.K,byw^, 



and in the second place we observe that 

-S-^ + Ll]W, 



dS 



dW 


/ dS 


dtj,a 


^dtj^a 


dW 


/ dS 


dtjM 


\dtj^b 



dW 

dtj,a 

dW 

dtj,b 



^dtj,a 

' dS 
^dtjb 



Li)W. 



Now, using the factorization problem we get 



dtj.a 

dS 



+ Li 



dS 

dtna 



s-\ 



dtj,b dtjb 

which, taking the + part (upper triangular) and the — part (strictly lower triangular) imply 



dS 

dt:,,a ' 



<Li)^, 



dS 
dtj„ 



s- 



(Lih, 



dS 



dt 



s- 



j,b 



dS 



dt 



s- 



j,b 



(124) 

(125) 
(126) 



{-Bl, + 5b,b'x')<'- (127) 



(128) 

(129) 
(130) 
(131) 



(132) 
(133) 



(134) 



so using (134) into (132| and (133) with the definitions (115) we obtain (124). The linear system (125) is obtained by 
inserting \112\ into (|124| . 
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To obtain the Lax equations (128) we take derivatives of (114 1 





r as 






dLa' _ 


r dS 


9tj,b 





-S , La' 



— [Bj,a, La'], 
= [BjM,La'], 



dtj.a 

dh' 



dt 



3,b 



dS 
dS 



S ^,Lb' 



dt 



3,b 



S ^,Lb' 



- [BjM,Lb'], 
~ [Bj^b, Lb']- 



Finally, (1291 are obtained as compatibility conditions for (124). 



□ 



All these equations provide us with different descriptions of a multi-component integrable hierarchy of the 2D Toda 
lattice hierarchy type that rules the flows of the multiple orthogonal polynomials with respect to deformed weights. 
This integrable hierarchy is the Toda type extension of the multi-component KP hierarchy considered in [9]. 

3.3 Darboux— Miwa discrete flows 

We complete the previously considered continuous flows with discrete flows, which we introduce through an iterated 
application of Darboux transformations [5] . 



Definition 20. Given sequences of complex numbers 

Aa := {Aa(?T-)}«ez C C, a = l,...,pi, Ab := {Ah(n)}„gz C C, 6=1,, 



,P2, 



(135) 



(where X is not intended to denote the complex conjugate of X) and two vectors, (si, . . . , Sp^) G and (si, . . . , Sp^) S 
Z^^, we construct the following semi-infinite matrices 



pi 



D. 



)b 



6=1 



n:"=i(AM 


- K{n)^i,a), 


Sa > 0, 


ni,a, 




Sa = 0, 


nL^='i(AM 


- Xa{-n)Wi^ay\ 


Sa < 0, 


nti(Ai, 


- Afc(n)n2,b), 


Sb > 0, 


n2,6, 




Sb = 0, 


nlf:'i(Ai6 


^Xb{-n)I{2,b)-\ 


Sb < 0, 



(136) 



where s := {sa, Sb}a=i,...,pi denotes the set of discrete times, in terms of which we define the deformed moment matrix 

fc=l,...,P2 

gis) ^ Do{s)gDo{sy\ (137) 
Proposition 24. The moment matrix g(s) has the same form as the moment matrix g but with new weights 



Wl,ais,x) = &aix,Sa)wi^aix), 



n^Lil^;- Aa(^^)), Sa>0, 

1, Sa = 0, 

'n?=i(^-AbH), sb>o, 



(138) 



W2,b{s,x) = %{x,Sb) W2,b{x), •= 1' 



Sb = 0, 



Thus, the proposed discrete evolution introduces new zeroes and poles in the weights at the points defined by 
sequences of A's. For example, in the a-th direction, the Sa flow in the positive direction, Sq — > Sa + 1, introduces a 
new zero at the point Xa{sa + 1), while if we move in the negative direction, Sa — > Sa — 1, it introduces a simple pole 
at Xa{sa — 1)- Let us stress that for the time being we have not ensured the reality and positiveness/negativeness of 
the evolved weights, this will be considered later on. 

3.3.1 Miwa transformations 

Here we show that the discrete flows just introduced can be reproduced with the aid of Miwa shifts in the continuous 
variables. 



Definition 21. We consider two types of Miwa transformations: 
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1. We introduce the following time shifts 



i T := \ij,a' T^a'.a^-^, ^j,b'\ j = l,2...., 

JZ-' a' = l,...,pi, 



(139) 



b =l,...,p2 



2. Dual time shifts are 



t^t±[z 1]^ := |tj,„,, t.i,b' ±db',b^\ 3=1,2, ■ 

^ a' = l,...,pi, 

b' = l,...,P2 

Proposition 25. The Miwa transformations produce the following effect on the weights 



(140) 



t^i,a'(a;,t T [z \,s) = (1 - 



Wi^a'{x,t± [z \,S) = Wi^a'{x,t,s), 



Wi^a'ix,t,s), W2,b'ix,t^[z ]a,s) ^ W2,b'{x,t,s), 



W2,b'ix,t± [z \,s) = 1 



Proof. When we consider what happens to the evohitionary factors under these shifts we find 
exp(^ij-a,a;^) ^ exp (t^-a' T (^a'.a^)) = (l^f)"^ ° '° exp ( ^ i^- 



(141) 

W2,b'{x,t,s). (142) 



(143) 



and therefore the weights transform according to 



Wi^a'{x,t,s) ^ [l - -j ' Wi^a'{.X,t,s), 



(144) 



which is hke the Darboux transformations considered previously. For the dual Miwa shifts we consider what happens 
to the evolutionary factors under these shifts 



exp I 



J,b' 

and the transformation for the weights is 



tj^b'X^j exp ( - X! y'3'''' ^ ^''''''^ 

J,b' 3,b' -'^ 



X 



±1 

1 ) exp 

z 



{-^^],b'X^ 
j,b' 



(145) 



W2,b'ix,t,s) 



X \ '^^b b' 

1 ' W2,a'{.X,t,s). 

Z / 



Thus, a comparison of (138), (141) and (142) leads to 



Proposition 26. Miwa transformations and discrete flows can he identified as follows 

CaWi^aix,t,Sa) 
CbW2,b{x,t,Sb) = < 



(wi^a{x,t-J2n=l[^a{n) ^ 




Sa > 0, 




-Aa(n))-i, 


Sa > 


< Wi^a{x,t), 




Sa = 0, Ca := < 






Sa =0 


[w,,a{x,t+j:^::},[Xai~n)- 




Sa < 0, 




-Aa(-n)), 


Sa < 0, 


'yj2,bix,t + j:lu [Mn)-'] 




Sb > 0, 


\UlUi- 




Sb>0 


< W2,b{x,t), 




Sb = 0, Cf, := < 






Sb^O 


[w2,b{x,t-j:^:^, [~Xb{-n)- 




Sb < 0, 


inli'i(- 


-Xb{-n)), 


Sb < 0. 



(146) 

□ 



(147) 



As a conclusion, the discrete flows and Miwa shifts in the continuous flows are the very same thing, and therefore 
we could work with continuous flows and Miwa transformations or with continuous/discrete flows. This discussion 
justifies the Miwa part in the name we gave to these discrete flows. 
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3.3.2 Bounded from below measures 

Of course, in order to preserve the link with multiple orthogonal polynomials, these discrete flows must preserve the 
reality, regularity and sign Constance of the weights, which generically is not the case. When the support of the 
weights is bounded from below, i.e. there are finite real numbers and iff,, such that supp(u;i.ad/i) C [Ka, oo) and 
supp(w2,f)d/i) C [^f,,oo), a possible solution is to place all the new zeroes and poles in the real line but outside the 
corresponding support, \a{n) < inf(supp(wi_ad/i)) and Xb{n) < inf(supp(w;2,bd/i)). A different approach, which will 
be considered in Appendix [B] is to arrange the zeroes in complex conjugate pairs. 

To analyze the consequence of the discrete flows on the integrable hierarchy we introduce two sets of shifts operators: 

Definition 22. 1. Let us consider the sets of shift operators {TalaLi o,i^d, {Tb}^'^^, where Ta stands for the shift 
Sa i-> Sa + 1 and Tj, stands for Sb i— >■ + 1. The rest of the variables {sa', Sb'} will remain constant. 



2. We introduce 



:=I-ni,„(I + A„(s„ + f)) + Ai,„, 

(148) 

96 I - n2,b(I + Af,(sfa + 1)) + Aj,. 



3. We also define the operators 



6a SqaS-^ = I - Ca{l + \a{Sa + 1)) + ia, 4 := SqbS-^ = I - Cb{l + M^b + 1)) + U, (149) 

Cq := S'lIi^a'S' ^, Cb '■= 5112,65 ^. 

Here the matrices 6a and 6b are called lattice resolvents. 
4-. Finally the semi-infinite wave matrices 

W := SDo, W := SDq. (150) 

Observe that 



iTaDo)D^' ^ Qa, Do\TaDo) 

Do{fbD^^) = qb, {fbDo)Do^ 



(151) 



When we assume that the semi-infinite matrices 6a and 6b are LU factorizable as in (33), i.e. all their principal minors 
do not vanish, we can write 

6a=6a^^6a,+ , 6b^6^^6b^+, (152) 



where 6a,- and 6b,- are lower matrices as is S* in (33), and 6a,+ and 6b,+ are upper matrices as 5 in (33). We now show 
that when the deformed moment matrix g{s) is factorizable, and therefore the multiple orthogonality makes sense, the 
following holds 

Proposition 27. // the deformed moment matrix g{s) is factorizable for all values of s then so is 6a and 6b with 

6a,+ ^{TaS)S~\ 6a,- ^ {TaS)S-\ ^ ^ 

_ _ _ , _ _ , (153) 

6b,+ = {TbS)S'\ 6b,+ ^ iTbS)S-\ 

Proof. When we apply the discrete shifts to the Gauss-Borel factorization problem g{s) = S~^{s)S{s) we get 

Ta{S~^)Ta{S) = Tag{s) = {TaDo)D^'g{s) = qag{s) ^ {{TaS)S-Y\TaS)S-^ = 6a, 

fb{S-^)Tb{S) = Tbg{s) - g{s)Do{TbD^') = g{s)qb ^ {{fbS)S-^)-\{fbS)S-^ = 6b, 

and the desired result follows. □ 

Therefore, we can consider the following 
Definition 23. The semi-infinite matrices uja and Cub are given by 

i^a-^ 5a,-6a^ 6a,+ , CJb 6b,- = 6b, + 6^^ , (154) 
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and show that 
Proposition 28. 



1. The following auxiliary linear systems 



TaW = UJaW, TaW = WaW, 

are satisfied. 

2. The Lax matrices fulfill the following relations 

3. The following discrete Zakharov-Shabat compatibility conditions hold 

{TaUJa')uJa ^ iTa'UJa)uJa', {TaOjb)uJa ^ {TbUJa)uib, {fbU)b')COb ^ {% U}b)t^b' ■ 

4-. When the discrete and continous flows are considered simultaneously, the following equations 



(155) 



TbBj^a — {djM^b)^b ^ + ^bBj^a^fj 

Tb'Bj,b = {dj^b^b')^b'^ +^b'Bj^b^i^,^, 



(156) 



(157) 



(158) 



are obtained. 
Proof. We compute 



T,W = iTaS){TaDo) = iTaS)S-'SqaS-'SDo = 6a,-SaW = Sa,+ W, 
TaW = {TaS)Do = {TaS)S-^SD„ - da, + W, 

TbW = iTbS)Do = {TbS)S-^SDo = 6b,-W, 

fbW = {TbS){TbDo) - {TbS)S-^Sq^^S-^SDo = Sb,+6b^W = 5b,-W, 



from where we deduce (155), which in turn imply (156) and (157) 



□ 

The simultaneous consideration of continuous and discrete flows leads to the replacement Wo — > WqDq and Wo — > 
Wq-Do, and the corresponding modification of the weight's flows is achieved by the multiplication of the continuous 



and discrete evolutionary factors, in this context we also have (158). These discrete flows could be understood as a 



sequence of Darboux transformations of LU and UL types in the terminology of 5J, which motivates the Darboux 
part in name we give to these discrete flows. In fact, we have that the lattice resolvents satisfy 



Sb = SZ^Si 



TaSa = UJaSaUJ^ ^ = 6a^ + S^^^_Sa, + Sa : 

fbh = WfcfJfcW^^ = 5b.-5b-^b,+5b^^ = 4,+<5f^i 



which amounts to the typical permutation of the LU factorization to the UL factorization. When there is only 
one component we have 5 ^ L + \ and (5 = A + L and the shift corresponds to the classical LU or UL Darboux 
transformations. 

If A^'a'^ (cc, s) , A^j^^ (x, s) are the multiple orthogonal polynomials and dual polynomials in the x variable correspond- 



ing to the discrete evolution of the weights ( 138 ) respectively we have the discrete version of Proposition 23 



Proposition 29. The wave and adjoint wave functions (150) are 

¥^\z, S) = 4^)(Z, S)3^a{z, Sa) (^D^'^C^, s) = (z, s)%{z, Sb)-\ 

and the the linear forms 



Q«(x,5)=£A«(x,s)z«i,a(x,s), 



(x, s) = £(^)('=)(x, S)w2.b{^, s), 



(159) 



(160) 



6=1 
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associated with the weights wi_a{x, s),W2^b{x, s), can be expressed as 



pi 



P2 



b=l 



in terms of which we have the equations 



Z ~ X 



-W2.b{x)dn{x), 



1 



-Wi^aix)dfl{x) 



(161) 



(162) 



Here the Cauchy transforms must be interpreted in exactly the same terms as in Proposition [s] Observe that ( |162[ ) 
do not correspond to the functions of the second kind 



Ci'\z,s) := 



QW(x,s) 



W2,bix, s)dfi{x), 



Cf\z,s) := 



Q«(x,s) 



Wi^aix,s)d^J.{x) 



(163) 



Notice also that from (149)-(152|, relations that hold true for any g and not only for the moment matrix, we get 



Lemma 6. We have that 



Co, = (^fc^o(A^)l"^|-"^'^ + l + ^,^i(A^)l"^|-"-'' + • • • + Qb.\n.\-n.,, + l, 

+ ■ ■ ■ + Pa,|ai|-ni,a + l, 



^T'j|n2|-n2.b + l _|_ ^(AT-j|n2|-ri2.i, 

Pa,o(A^)l"^l-"i-+i + Paa(A^)l"il-"^- 



— 7,, A|n2|-n2,b + l I A|n2|-n2,b i i „ 

— PbfilV -\- Pb^llV + +P&,|n2|- 



■n2,!, + l' 



for some diagonal semi-infinite matrices 



with 



a J ^ diag(a;Qj(0), Waj(l), . . . )■ 
u)b,j = diag(wbj(0),W6,-,(l), . . . ), 

Pa J = diag(pa,i(0), Paj(l), ■■■), 

Pbj = diag(pf,j(0),p6j(l), . . . ), 



Paj(fc) := '^aj(fc -\ni\+ ni^a - 1 +i), 

P6j(fc) := '^6j(fc + |?l2| - »^2,fc + 1 - j)- 

that with 
Definition 24. M^e define 

7a,a'is,x) := (1 - (5a.a'(l + Aa(.Sa + 1) - 2:)), 

76,6' (s, a;) (1 - (56,6' (1 + A6(s6 + 1) - a;)), 

leads to 

Proposition 30. The following equations 

A« = -.vo(fc)Ar''"'^"^-^'"'^ + • • • + ^.',i«2i-n2,b'+i (fc)^i'^- 

p,',o(fc)(T„'4'^i"^i+"^-'~'^-...-^ . ^fi^v^.m^- m 



+ • • • + PaMft,|-„,,„'+l(fc)(Ta'^r) = < 



(p6',o(fc)(Tl 



6' A 



(fc+|n2|-n2,b'+l) 



) + •••+ PbMn2|-n2.b'+l(fc)(r&'4") = A"-; 



are fulfilled. 



Proof. For ( 168 ) recall the discrete auxiliary systems for W , while for ( 169 ) just consider that 



1\T 



(164) 



(165) 



(166) 



(167) 



(168) 
(169) 



□ 



Notice that relations ( 168 ) and ( 169 ) are among multiple orthogonal polynomials in the same ladder but with 
different weights, they link the polynomials for the weights Wi,^, W2,b with those with Ta'Wi a, Ta'W2,b or T^iWi^ai TyW2fi. 
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3.4 Symmetries, recursion relations and string equations 



We now return to the discussion of the symmetry of the moment matrix that we started in ^2.6 but with evolved 
weights and the use of Lax matrices. The first observation is the following 



Proposition 31. The j-th power of the evolved Jacobi type matrix introduced in ^2.6\ is related with Lax matrices 
through what we call a string equation: 

Pl P2 

^Y.La^Y.U>^ J = 1,2,..., (170) 

a=l b=l 

and the multiple orthogonal polynomials are eigen-vectors: 

J^K'^x^K', (J^y^b. =xWfc,, (171) 

for a' = 1, . . . , Pl and 6' = 1, . . . , p2 • 



Proof. Using ( 68 ) it can be proven by induction on j that for any j > 1 the following equation holds 

Al,gU2,t = Ili.ag{Alf,y, (172) 

so that 

Lia = Call (173) 



Summing over a,b we deduce (1701. Moreover (1711 is obtained as follows 



Pl 



.P^a' = SY^KaS^'SxiX = X^^a', (174) 
a=l 

P2 

{J^Vda. = iS-YY.^iaS^(^'yX2,b' = XW,,. (175) 
b=l 

□ 



We are ready to show that the symmetry ( 68 ) induces a corresponding invariance on Lax matrices and multiple 
orthogonal polynomials 

Proposition 32. The following relations hold for j ~ 1,2, . . . 

a=l f)=l J'" a=l -I'"- b=l 

Pl Q P2 C) P^ ^ P'^ Fl 

a— 1 b—1 a— 1 b—1 

Proof. See Appendix |A] □ 
3.5 Bilinear equations and r-functions 

The proof of the bilinear identity needs three lemmas. For the first one, let Wfi^^fi^,Wfi^^fi^ be the wave matrices 

associated with the moment matrix gni,n2j so that, Wa^fi^gft^fi^ = Wft^^^. Then, we have 

Lemma 7. The wave matrices associated with different compositions and times satisfy 

Wn,,nAt.^)^\,nWn',,n'St'^sT^ = W^3i (^ , s)7r J, (t' , s')- \ (178) 

Proof. We consider simultaneously the following equations 

Wft^ M2 {t, s)g = Wn^ ,n2{t,s), 

where g = gnt,n2, and we get 

Wfli,fj2(i,s)"^W3i,,-j,(i,s) ttT, -^W^f;'^,H^(<',s')"^PFa;,s^(i',s')7rf;^,n2 =9, 
and the result becomes evident. □ 
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For the second one, let (•)-! denote the coefficient in z ^ in the Laurent expansion around z ^ oo (place where 
the Cauchy transforms make sense). 

Lemma 8. For the vectors Xa the following formulae hold 

(Exa(x:n_i = (Ex:xi)_i = i, 

a— 1 a— 1 

and therefore 

Lemma 9. For any couple of semi-infinite matrices U and V we have 



pi 



UV={J2iUxi,a){V^Xla) 
a—1 



(179) 
(180) 



b=l 



Proof. It follows easily from Lemma [8} 



{j2^Uxi,a){V^Xl,aV)_^^u[j2xiAxlaf)_V^UV, 
a=l a—1 

{T.(^xh)iv^X2,bV)_^ = u{J2xhxlb)_v - uv. 



6=1 



6=1 



□ 



We have the following 
Theorem 7. 1. The wave functions and their companions satisfy 

Pl r P2 



a=l"'°o 6=1"'°° 

2. Multiple orthogonal polynomials, their duals and the corresponding second kind functions are linked by 

P^ I- P2 „ 

E f 4':U.a(^,^:^)4i..i.a(^-^'>^')i?a(.)dz = f ^^Sk .6 ^ : ^)4! .6 ^' > ^')^6 (.)dz, (181) 



whe 



Ea i^a^a){z,t,s){{r^a^a){z,t',s'))-\ 



Eb {ib^b)iz,t,s){{ib^b)iz,t',s')y 



Proof 1. If we set in ( [T79| ) U = Wft,,n2i't, s) and V = 7rT^.fj^VKff^,fli(t', s')"^ and in ( |T80| ) wc put U = Wft,,n2it,s) 
and V = ■^^.n^^n'i^n'ii^' ^ ^')~^ attending to ( |178[ ), recalling that ^'Si,S2,a = Wii,,n2Xni,a, *ni,n2,b = ^ffi,n2Xs2,6 
and observing that = (VF-^^'«^)^^s;,a,XHi,a and *t,.Ki,6 = ^S2X«2,6 we get the stated 

bilinear equation for the wave functions]^ 

2. We can write 



which strongly suggests to consider in (179) 



^The reader familiarized with Toda bilinear equations should notice that in the right hand term we are working at 2: = oo instead of, 
as customary, at 2; = 0; the reason is that for the definition of J24 we have used X2 instead of X21 order to get polynomials in 2, while 
normally one gets polynomials in . See footnote 2. Moreover, due to the redefinition of the wave functions there is no factor 
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Analogously 



suggest to set in ( 180 ) 



The application of ([l79|, p^ a nd ( [T78| gives the alternative bilinear relations (flSll) where we have used the 
evolved Cauchy transforms ( 123 1 and introduce the evolutionary factors 



Ea ■■= {'ga&a){z.t,s){{<Sa&a){z,t' ,s'))-\ 



□ 



The factors involved in this definition were introduced in (109) and (138), so that we assume the discrete flows 



within the bounded from below support scenario, while if we consider the two-step discrete flows the replacement of 



the ^-factors by the factors ( 220 ) is required 



It can be shown that for certain weights, for which we can use the Fubini and Cauchy theorems, and when one 
only considers a finite number of continuous flows that the r.h.s and l.hs. in this bilinear relations are proportional to 
Jg^Q^^^ fi^ix,t)Q'"S! ^, {x,t')d^{x). This is a direct consequence of 

Proposition 33. We have the following identity 



Proof. See Appendix [X] 



(182) 
(183) 

□ 



Now, we will perform a full characterization of the r-functions associated with the multiple orthogonal polynomials 
defined in this paper. 



Definition 25. Let us define the following matrices 



( 90,0 90,1 
91,0 91,1 



9+a 



9o,i-i 90,1 \ 
51, (-1 9i,i 



91-1,0 91-1,1 ••• 91-1,1-1 9i-i,i 

\~9l + a,0 ~9l + a,l 9l + a,l-l 9l + a,l J 



I 90,0 90,1 
91,0 91,1 



y+b ■ ^ 



.90, /-I 
9i,i-i 



9i-i,o 9i-i,i ■ ■ ■ 9i-i,i-i 
\ 9i,o 91,1 ■■■ 9i,i-i 



9o,i+b \ 
9i,T+b 



9i-i,i+b 
9i,i+t ) 



(184) 



The matrix 5!^^^' is obtained from .g^'+^l replacing the last row (operation denoted by a dashed line) by 

(.9J + a,0,5i + a,l' ■ • • T9l + a,l-1^9l + a,l)^ 

and 5+t^'is obtained from eplacing the last column by {90,1+^^ 9i,i+b' ■ ■ ■ ' 9i-i,i+b' 9i,i+b^^ ■ clear that if 

ai(0 = a then .g![+^l = .gl'+^l and if 02(1) = b then gl'l^' = g^^+^l 



The minors of the these matrices ( |184| will be denoted as m|'+^1 = MH'' for gl'+^l, M^^^^ for g^l"' and MH-'.^ 
for 5!^^^'. Now we introduce the following determinants that are cofactors of the previously defined matrices 

Definition 26. The r-functions are defined as follows 



r['+l] 



M+1] m['+i1. 



.(0 



' , ■■= (-1) 
a,— a' \ / 



'+b,-b' 



~--a -b - (-1)'-"+'-"^,''+'' 



l-a,l- 



(185) 
(186) 



Moreover, 
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1. If ai{l) — a then we denote ri'^, := r^'^ and r/''' 



.(0 



b ■— '-a,-b- 

2. If a2{l) = b then we denote ri'^ := tI'^ _^ and f^l, := f^l _j^,. 

3. We also introduce r^'^ = f := det gl'l and 

'+a ■— "-^^^ y+a ' 

If ai(0 = a then t|'+^^ = t('+i\ and if a2{l) = b then ff+^^ = r 



:=det^!^^l 



Given a perfect combination (/k, wi, and the corresponding set of multiple orthogonal polynomials {^[£7i:,72].a}aLij 
with degree vectors such that It'll = |i72| + 1, there exists a (ni,r?2) ladder and an integer / with |i7i| = I + 1 and 
|j72| = I such that the polynomials {>la'}aLi coincide with {^[i7i;i72],a}aLi- The final result does not depend upon the 
particular {Hi, 712) ladder we choose to get up to the given degrees in the ladder; however, the r-functions do indeed 
depend on the ladder chosen through a global sign. A simple sign-fixing rule is to choose the ladder ni = ui and 
n2 = + . We define 



= \Vi\ - 1 = |z72|, 



and we deduce 



Proposition 34. Given degree vectors [vi, 1/2) such that \vi \ = \i'2\ + 1, « composition with n\ = Vi and n2 = 1^2 + ^P2 
and ^ = It'll — 1 = |t72|, we have the following identities 

T+l_a' = a')rp^-ei,„,+ei,,;i72]' ^+5,-6' = ^2,2(&, 6')Tp^;i72-e2,t'+e2.!,] ' 

(0 _ -(0 _ (K \ 

T-b-a — '^-a-b — ^2,1 (,0, ajT^^+^j ^^ .g'l „;i72+e2, 



where 



|-e2,P2~62,6j' 



In particular 



£i,i(a,a') 


:= (- 




a' 


< a 


ei,i(a, a') 


:=(- 


_l)E°=l "l.j+IZ^Ll "Ij+Sa'.pi ^ 


a' 


> a 


£2,2(6,6') 


{' 


_t)E5=i ''2,3+EjLi "2,3-1^ 


b' 


< b 


£2,2(6,6') 


■■= (' 


_1)E5 = i ''2.j+E5Li '^2,j^ 


y 


> 6 


£2,1(6, a) 


:= (- 


_]^^Ej = l ''2,j+E°=l 1'1,J+5!),P2 ^ 






£i,i(a,a) 


:= 1 


= £2,2(6, 6). 







T^a = £l,l(Pl>a)7'Pi+ei,„-ei,„;i72]> 
£l,l(a,Pl)T'[i7i+ei,a;i72+e2,p2]' 



T^-l = £2,2(p2,6)Tp^;i72+e2, 



P2 ~^2,b\ ' 



'+a 



r+b = £2,2(6,P2)r[i7i+ri,p^;i72+e2,6]- 



We now proceed to give the r-function representation of multiple orthogonal polynomials, their duals, second 
kind functions and bilinear equations. The r-functions allow for compact expressions for the multiple orthogonal 
polynomials: 

Proposition 35. The mixed multiple orthogonal polynomials Aa\ a ^^'^ ^-la ^'^^^ following r-function 
representation 



•^1.0.(1)- 



T(0(f) 



4(n,a') 

^Pi(i-l)+ei,„,;i72(i-l)],a 



[i'i(l);f'2(l)-e2,b],a 



r('+i)(t) 



I > 1, 
I > 1, 
I > 1. 



(187) 
(188) 
(189) 
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The dual polynomials A^''', A^j^^j^, and ^ have the following T-function representation 



Proof. See Appendix [X] 

Observe that in the simple ladder defined above {i^i, V2 + (i2,p-,) with Z = |i72| = l^'il ^ 1 we have 

vi{l) = vi, 1^2(1 -i) = V2, 

vi{l -I) = vi~ ei,pi, V2[l) ^V2+ e2,p2. 

From Proposition [35] we get 

^(II,Pi) /X f \ ^ i'''pl+ei,p -ei,„a;i72](^ ^ "^]a) 

A....l.W = ^M(Pi,a)z- ^—^^ , 

..;.2+e2,.2-e.,.]-^ 2,U rp,+e,.,,;.-2+e-2.,,]W 



^P2+e2..';^^i-e.,.J,^ - ^2,2l^>,(^ ^^^^^ ^^^^^ 



A-' 



(I,a) 



e-l,a],f> -^2,1(6, aj^ — ^ — -T-. • 

' |i^i+ei,pj ;!^2+e2,P2l ' 



We now present the r-representation of the Cauchy transforms of the linear forms. 
Proposition 36. The Cauchy transforms have the following T-function representation 

r('+i)(i) 

rW _ -.-'^2.!,(;-i)-i '^+b ~ ^^b) 

r(0(t) 

Proof. See Appendix [X] 

We have the representation 

db'^^l - . , (z) ^ £1 i(a,pi)z~-^^--i+^"-^i ^P.+e-:..;.-2+e-2,P2l(^ + [^'']°) 

[.2+e2,P2;-i-ei,pJ,a^ ^ 1^1^ '^1^ +e,,p, ;.~2+e2.P2 1 W ' 

dl'-l^] , (^) ^ S2 2{b,P2)z--^'^-' ^'-"^^'^^■-^^--+-----'/^^ ~ [^''^"^ . 

Finally, we consider the r-function representation of the bilinear equation 
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Proposition 37. The r functions fulfill the following bilinear relation 



E 

b=l 



X.6(0-'-2,t,(fe-l)-2-('c+l) 



- [--']t)rl!U-t(:^ + [z-\)E,iz)dz. (197) 



Proof Just consider (181) together with (187), (190), (195) and (196). 



□ 



This bihnear relation can also be written as follows 



^ J Z — QQ 

z 



^'^2,b+^b,p2-'^2,b-2 



(198) 



That with the identification m* = vi+ ei ^^, n* ^ 1/2, rn = iy[ and ri = + 6*2.p2i to signs, is the bilinear relation 
(41) in [g. 
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Appendices 

A Proofs 

Proof Proposition [7] The orthogonality relations can be recast into two alternative forms 



1,1 



( .90,0 30,1 
.91,0 ffi,i 



90.1-1 \ 
91.1-1 



= - {91,0 9i,i ■ ■ ■ 9i,i-i) 



Si.i-i s, 



\5i-i,o ffi-i.i ••• 91-1,1-1/ 

/3o,o 5o,i ••• 5o,i\ 
51,0 51,1 ••• 9i,i 



i.i I 



(0 ••• s 



1,1 ) 



\9i,o 9i,i ■■■ 9i,iJ 



I + 1 components 



From ( 43 1 we get 



(fc) 



A:=0 



= ^['^-{m,o 91,1 ■■■ ffu-i)(5['i)-^ef' 

= 5u(0 ••• I) {gV+'^r'^f+'l 



I > 0. 



use ( 199 ) 



use ( 200 ) 



(199) 



(200) 



(201) 
(202) 
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Cramer's method solves ( 199 1 as follows 



det 



where Mf^^ is the (i,j)-minor of the truncated moment matrix (7W defined in (21). Therefore, 



i=0 





' 50.0 5o,i • 


■ 90,1-1 


^(0) 


jTj- det 

det gV\ 


51.0 514 • 


■ 91J-1 


Ai) 
?i 


3;-i,o • 


■ 9i~i,i-i 


Ai-i) 




V 91.S) 9L1 ■ 


■ 9LI-1 


^(0 



I > 1. 



The orthogonality relations for the dual system can be written also in two alternative forms 



/ 5o.o 5o,i 
51,0 51,1 



9od-i \ / -^o,; \ 



91J-1 



\3;-i,o 91-1,1 ■■■ gi-1,1-1/ \Si-i,i/ 



(90,0 9o,i •■• 9o,i\ 
9i,o 9i,i •■• 51.; 



/^o,A 



\3;,o •■• 9i,iJ \SlJ 











/ 5o,; \ 
9i,i 

\9i-i,i) 



I > 1, 



I > 0. 



As before, (43) leads to the following expressions for the dual linear forms 



)(0 



2 



k=0 



( 9o,i \ 











51,/ 



\9i-i,i/ 



use (2041 



use (2051 



From ( 205 ) we obtain 



(_l)i+JM/y 
det(<7['+il) 



j=0,...,l, 
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and consequently 



J=0 



2 



3=0 



det 



det 



/ 9o,o 50,1 • • • gaj-i 

91,0 giA ■■■ gid-i 



91-1,0 9i-i,i • • ■ 91-1,1-1 



S2 



.90,; \ 



I > 0. 



Proof Proposition [9] We have 



p(0 ^ 

detgW 



fe+iA,f(' + l) 
k,l 



E 



fc=0 



fe2=f2,6('-l) 



a;'''^''-'^«i,ai(fe)(2^)^^2,6(a;)x'=MAi(a;), 



which according to ( 20 ) recasts into 



det 



fe=0 





^ 90,0 


90,1 


■ ■ 9o,i-i 


f(0 

^ 0,6 


FTT det 

det 


91,0 


91,1 


■ ■ 9i,i-i 


f(0 

^ 0,6 


91-1,0 


91-1,1 


■ • 9i-i,i-i 


f(0 

^ l~l,b 




y 91,0 


91,1 


■ ■ 9i,i-i 


f(0 



We also obtain 



\l+kj.M+l) 
fc=0 ki=vi,a{l-l) 



-ki-1 



which can be written as ( 20 1 



^'^ - detg['+il ^^^'■'^ ' 



k=0 





/ ffo,o 


9o,i ■ 


• 5o,/-i 


.90, i 


tTttt det 

det 


9i,o 


9i,i ■ 






9i-i,o 


9i-i,i ■ 


• 9i-i,i-i 


91-1,1 








-^/-l,a 





Proof Proposition 15 From (53) and (54 1 we deduce that 



Uz))^<,(^') = (xlaW)^Xl,a'(^') 



Z — z' 
5k 



{^,{Z)Y ^'Az') - (X16(^))^X2,6'(^') - ^ 

Z — ; 

[^a{z))''%{z') = (x^,a(^))^5X;,6(^')- 



; > 1. 



I > 1. 



z' < z 



z' < z , 



(209) 



(210) 
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The two first relations imply the corresponding equations in the Proposition. For the third we observe that from ( 10 ) 
we get 



= I (Xi,q(^)) XiA^)iX2,b{x)) X2,6'(^')'^i.a(2;)w'2,6(a;)d/i(a;) 
wiA^)w2,b{x)dn{x) 

wi,a(a;)w2,b(a;)d^(a;). 



{z — x){z' — x) 

I r / 1 1 



z — z' \z — X z' — X 



Proof Proposition 20 Using Definition |l7| for the linear forms Q^^^ and multiplying by (dk^))^ have 
integrating both sides we get 



Q^iW(4"(^))^dM(x)= / £,^|^-\x){^'i{x)VA^i{x) 



= 0, 

that written componentwise gives the following orthogonality relations 

/ Q^i(x)w;2.a.(fc)(a;)a;'^*'MM(a^), fc = 0, . . . , / - 1, 

or equivalently 

Ql\{x)w2,b{x)x^Mx) = Q, < fc < j.2.f,(^ - 1) - 1, 6=1,...,P2. 



Notice that, ^ is monic and degA^|j^(x) — ki{l+a) but A^j^ ^, with a ^ a' satisfy degA'j'lj^, < ki{{l — 

',(0 



This means that the set of polynomials A^a'i^) have degrees determined by — 1) + ei^a and a normalization 

M) _ ^(II,a) 

!+a — V[£7i(;_i)+gi^;,72(;_i)]- 



with respect to the a-th component of type II; i.e, = Q^V'"^ 



In a similar way, the associated linear forms (5^)j(x) solve a mixed multiple orthogonal problem that can be obtained 
as follows. From Definition [l7| and multiplying by (C2'(^))^ '^^ S^t 

integrating both sides 

n{x){i^'\x))A^,{x) = eJ^g^^+'Y' I il'-'\xM^'\x)YMx) = ^, 
and written componentwise 

that is equivalent to 

A^{x)w2,b{x)x''Mx)^5^j_^, < fc < ^2,&(0 - 1, 6=1,...,P2. 
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Hence, the set ^^j^ ^, is a type I normalized to the 6-th component solution for a mixed multiple orthogonality problem; 
,(0 



the degrees satisfy degA_j,^, < l-a'- Moreover, the fact that the last orthogonality condition in the 6-th component 
is missing gives the identification Q*''', — q\1:''),, ^ ^ 
Using Definition [it] and multiplying by (x) we have 



and integrating both sides 

I^Q%{x)U^f^dfi{x) = ( I ^f\x)^i+'\x)df,{x)- I ^f\x){^^\x)yd^,{x){g^'^)'^ 



^f\x)^^^^'\x)d^,{x) - {gm){g^'^) 



9i.i+i 



5i.r, 



( 9oJ+, \ 
9iJ+b 



\9i~id+b/ \9i-i,i^J 
= 0, 

that componentwise leads to the following orthogonality relations 

/ Qllix)wl,a,ik)ix)x''''^''^d^lix)^0, fc = 0,...,?-l, 
Jr 

or alternatively 



s(0 



ix)wi^aix)x''d^{x) = 0, 



< fc < i^i^l - 1) - 1. 



a = 1, . . . ,_pi. 



Notice that, ^^),^f, is monic and degA'df, = k2{l+b) but A^),;,, with 6 7^ 6' satisfy degA^),,,, < ki{{l - This 
means that the polynomials A^)^ j^, have degrees determined by V2{1 — 1) + 6*2.6 and a normalization with respect to 
the 6-th component of type II; i.e, Q+I = ^fv.Ai~i)+e^^b-Mi~i)y 

Finally, we obtain the orthogonality relations for the linear forms Q^^^{x). From the definition we get 

^i^'\x)Q^i{x) = cr^'(x)(er^i(-))"(5['+^i)-^e._„, 

and integrating both sides 



i^''i^x){i+'\x)Yd^^{x)W^+'y)-'el_^ = ei 



^^i+'\x)Q^!^,{x)d^i{x) -- 
and componentwise that means 

/ Qf!i(a:)a:'=iWu;i,,,(fc)(a:)d/i(a;)=4,i_„, fc = 0,---,;, 

or equivalently 



Q''l'^{x)x^wi^a{x)d^ji{x) = 6k,l_^, 



0<k< Vi,a{l) - 1, 



a = 1 



so the set A^^^^ is a type I normalized to the a-th component solution for a mixed multiple orthogonality problem. 
The degrees satisfy degA^_^\ ^, < l^i,'', and therefore we conclude that q''], = Q''^^^^y-.^^^^_g^ y 
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Proof of Proposition 32 Taking + and — parts in (170) we obtain 

Pl Pl Pl Pl P2 



a=l a=l 



a=l 



f)=l 



b=l 



b=l 



Using Lax equations and observing that LaLa' — La' La and Li,Lh' — Lb' Lb we have the following symmetries for 
the Lax operators 



d ^ d 



( E + E 5^)^''' - [ E B:,a + E La'] = EK, ia'] = 0, 

a— 1 /i — 1 



Pl 



6=1 
P2 



a=l 
P2 



; E + E a^)^"' - [ E + E B.b. Lb] ^ Et^i = 0- 

a=l ^'"^ f,= l o=l 6=1 6=1 



From (1271 we conclude that the multiple orthogonal polynomials and their duals are also invariant 

Pl Q f) 

( E ^ + E 3r^<' - (E + E B.-^ - = {J' - ^'X- = 

a— 1 "'^ 6^1 a— 1 6—1 



a— 1 b— 1 



6=1 



Proof of Proposition |33| We just follow the following chain of identities 

= Sn,fiAt,s)WoMAt.s){ j 6u(x)eT^(x)dAi(x))(W^o,«i(t',s'))~i(5B;,^Si(i',s'))~' 
= Sn^,n2{i,s)(^ j ^ni{x,t,s)^ff^{x,t',s')dfi{x)^{Sn'^^n'^{t',s'))~^ 



where (x, s) and fg' (x, t' , s') represent the vectors of weighted monomials but with evolved weights. 



Proof of Proposition 35 To find the r-function of the multiple orthogonal representation we first need two lemmas 
Lemma 10. Let R^^'> be the j-th row ofr'^'^t) and Ri^^ the j-th row of T^'-^t - [z~%), then 

(211) 



where j' = j + 1 if ri{j) < ni.a — 1, but j' = j + {\ni \ — ni_a) + I if = — 1. This is also valid for T_a 
and for r'fl . 



(0 Ai) 



Let now be C^^^ the j-th column of r^^^ and Cz the j-th column ofT^^\t^ [z ^]^) , then 



(212) 



where j' = j + 1 if r2{l) < 7i2,6 ~ 1 but j' ~ j + (1^21 — ^2,6) + 1 if f2{j) = ?^2,6 — 1- This is also valid for 
and for r^^ 



Proof. It follows directly from ( 141 ) and ( 142 ) 



□ 
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Let us recall the skew multi-linear character of determinants and the consequent formulation in terms of wedge 
products of covectors. Observe that 



Lemma 11. Given a set of coveetors {ri, . . . ,r„} it ean be shown that 

n n+l 



f\ {zr, - r,+i) = ^(-l)"+i-JzJ-i n A r2 A • • • A f, A • • • A r„+i, 



(213) 



n+l ■ 



□ 



where the notation fj means that we have erased the covector rj in the wedge product ri A • • • A r 

Proof. It can be done directly by induction. 

The proof of Proposition |35] relies on Lemma [TO] Lemma [TT] Corollary [T] and Proposition [22j First let's focus on 
( |187| ); it is clear that Vi'^(t - [z-^]a) expands in z according to ( |211[ ) for r^l and to ( |213| . Now n = ki{l^a) 

and the covectors rj should be taken equal to those rows i?'-'' with ai{j) — a. Observe that there are only ki{l_a){= 
^i,a{l) ~ 1) rows which are non-trivially transformed. In this form we get the identification of (49 1 with ( |187[ ), where 
the terms corresponding to the wedge with one covector deleted corresponds to the minors AlJf^^K Now, looking to 

(1881 and (1891 we expand again in z and use the same technique based on (211) for r^'^ and t1'^ _^ and (213 1. 
These allow to link pM to and (fT06| to pM. 



To prove (1901 we proceed similarly. Looking at (212) for rfj^ observe that there are only k2{l-b)i— V2fi{l) — 1) 
columns which are non-trivially transformed. Now, recalling (50) and using (213) but with Vj being the columns C'^^\ 
such that a2(j) = 6, and n = k2{l-b), we get the desired result. Finally for (1911 and (192 1 we expand around z to 
see the equivalence between ( 107[ ) and (1911 and the equivalence between (1051 and (192 1. 



Proof of Proposition |36| We need the following two lemmas: 

Lemma 12. Let i?^-'' be the j-th row o/ gl^^^' and R^P the j-th row of g^^'^^^t + [z~^]a), we get 

oo 
J' = l 

Let C^^^ be the j-th column of gi^^^ and Cz^^ the j-th column of c}j^^\t — [z^^]ff), then (142) gives 
Proof. For the first equality insert the expansion 



(214) 



i'=i 



(215) 



1 - 



OO U 



E 

fc=0 



into ( 142 1 . The other equation is proven similarly. 
Lemma 13. The following identity 

n oo 



ri A • • • A r„ 



3 = 1 j=0 



oo 
i=0 



holds. 

Proof. Use induction in n. 



□ 



(216) 



□ 



Finally Proposition 36 is proven using (210) and (209). 
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B Discrete flows associated with binary Darboux transformations 



When the supports of the measures are not bounded from below, (in which case the new "weights" ( 138 1 do not have 
in general a definite sign and therefore should not be considered as such), there is an alternative form of constructing 
discrete flows which preserve the positiveness/negativeness of the measures. The construction is based in the previous 
one, but now the shift is the composition of two consecutive shifts associated with the pair Aa(n) and Xa{n + 1), being 
complex numbers conjugate to each other; i.e., we consider 

Definition 27. We define a deformed moment matrix 

g{s) = D',{s)gms))-\ (217) 

with 

pi 

a=l 

[n:Li(|Aa(n)pni,,-2Re(A„(n))Ai,,+A?^J, > 0, (218) 

^0,a ■= \ ni,a, Sa = 0, 

lnL='i (|Aa(-n)|2ni,, -2Re{Xai-n))Aua + A?,J"\ Sa < 0, 



P2 

(5(,)-i:=E((I)„(.)')-i),, 

6=1 

f mUi\Mnrn,,b - 2Re(A,(n))Aj J + (Al^f)) 5, > 0, (219) 

((5o(s)')"')„:= <n2,b, Sb = 0, 

mlU {\M-n)\'ll,^b - 2Re(A,(-n))Aj,) + {Al,ry\ h < 0. 

Proposition 38. The previously defined deformed moment matrix corresponds to a moment matrix with the following 
positive/negative evolved weights 



W1,q(s,x) = ^'^{x,Sa)wi,a{x), 



W2,b{s,x) ^ ^l,{x,Sb) ^W2,b{x), {%) ^ ■■ = 







-A„(n)|2, 


Sa > 0, 








Sa = 0, 




\x 


-\a{-nT^ 


Sa < 0, 




\x 




Sb > 0, 


1, 






Sb - 0, 




\x 




Sb < 0. 



(220) 



Proceeding as in the previous case 
Definition 28. We introduce 

q'a:=l- Til.ail - \XaiSa + 1)|') - 2Re(Aa(Sa + l))Ai,, + A?,„ 

:= I - n2,5(I - \Msb + 1)1') - 2Rc(Af,(,Sfc + l))(Aj,) + (Aj,)^, 

and 

S'^ ■.= l~Ca(I~\XaiSa + l)\^)-2Rc{Xa{Sa + l))La+Ll 

S',-=i~ Cbil ~ \Xbisb + 1)1') - 2Re(Ab(sb + l))Lb + LI 
The wave and adjoint wave functions now have the form 

¥,^Hz,s)^Ai'Hz,sWz,Sa), {%)^'Hz,s)^A\;'\z,smz,Sb)-\ (223) 

and the expressions ( 161 )- p^2| still hold. 



(221) 



(222) 
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If we introduce lu'^ and a;[, as in (154) but replacing 6 by S' , the equations (155l-(158) hold true by replacement of 
u! by uj' . Now, the form u' differs from (164) as now we have 



^a-^a.O^ H ^ ^a,2(|«i|-ni,, + l)' 

4 = 4,o(A^)'""''""'''+'^ + -- 



With the definition of 



'"b,2(|n2|-n2.b + l)- 



j',y{s,x) := (1 - Sa,a'il ~\X- XaiSa + 1)1'), 

7ib,{s,x) := (1 - SbM'il Msb + 1)1'), 



(224) 



(225) 



we have that 



Proposition 39. The present setting ( 168 ) and ( 169 ) are replaced by 



k) 



(fc-2(|«2|-n2,,+l)) 



(fe-2(|ni|-ni „,+!)) 



) + ---+Pl'.2(|ff,|-„,,„,+l)(fc)(ra'4'^) 



A 



(fe) 

b ' 



,A 



(fc+2(|n2|-n2 .'+1) 



) + •••+ P6',2(|«2|-n2..'+l)(^^'4'^) )7M' = 4' 



(226) 
(227) 
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